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Abstract. We construct a theory of distributions in the setting of analysis on post-critically 
finite self-similar fractals, and on fractafolds and products based on such fractals. The re- 
sults include basic properties of test functions and distributions, a structure theorem show- 
ing that distributions are locally-finite sums of powers of the Laplacian applied to con- 
tinuous functions, and an analysis of the distributions with point support. Possible future 
applications to the study of hypoelliptic partial differential operators are suggested. 



1. Introduction 

The prevalence of fractal-like objects in nature has led both physicists and mathemati- 
cians to study dynamic processes on fractals. One rigorous way to do this on post-critically 
finite (p.c.f.) fractals is by studying differential equations in the natural analytic structure. 
A brief description of this analytic structure will appear in Section[2]below, but we empha- 
size that it is intrinsic to the fractal, and is not necessarily related to the analysis on a space 
in which the fractal may be embedded. For example, the familiar Sierpinski gasket fractal 
SG is often visualized as a subset of R 2 , but restricting a smooth function on R 2 to SG does 
not give a smooth function on the fractal [3 1. Similarly, we should not expect the solutions 
of differential equations on fractals to be quite like the solutions of their Euclidean ana- 
logues; for example, many fractals have Laplacian eigenfunctions that vanish identically 
on large open sets [2 |, whereas eigenfunctions of the Euclidean Laplacian are analytic. 

Perhaps the most important tools for studying differential equations in the Euclidean 
context are Fourier analysis and the theory of distributions. Since the theory of analysis 
on fractals relies on first constructing a Laplacian operator A, it is unsurprising that quite a 
lot is known about the fractal analogue of Fourier analysis. In some interesting cases the 
spectrum and eigenfunctions of the Laplacian are known explicitly, and many results about 
Laplacian eigenfunctions have also been derived by using probability theory to study the 
heat diffusion on fractals. Fourier-type techniques have also been applied to treat smooth- 
ness in the fractal setting: analogues of the Sobolev, Holder-Zygmund and Besov spaces 
that are so important in Euclidean analysis of differential equations were introduced and 
investigated in [33 1 . Analogues of other basic objects in Euclidean analysis are studied 
in fl9l l4l. By contrast there has not previously been a theory of distributions on fractals, 
and it is the purpose of the present work to provide one. 

It is relatively elementary to define distributions on fractals; as usual they are dual to the 
space of smooth functions with compact support, where a function u is said to be smooth 
if A k u is continuous for all k e N. The main theorems about distributions are then really 
theorems about smooth functions, and the key to proving many of them is knowing how 
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to smoothly partition smooth functions. Partitions of unity are used to achieve this in the 
Euclidean setting, but are not useful on fractals because products of smooth functions are 
not smooth [3|. (This latter fact also implies that products of functions and distributions 
are not distributions, so the distributions are not a module over the smooth functions.) 
Instead we rely on a partitioning theorem for smooth functions proved in [27], see Theo- 
rem l2.7l below. Using this partitioning result we are able to prove analogues of the standard 
structure theorems describing distributions as derivatives of continuous functions (Theo- 
rem I53, and identifying the positive distributions as positive measures (Theorem 15.61 ). 
We can also characterize the distributions of point support as finite linear combinations 
of certain "derivatives" of Dirac masses that can be explicitly described (Corollaries 16.61 
and 16. 13b , provided we make certain assumptions about the point in question. These as- 
sumptions are needed in order to understand the local behavior of smooth functions at the 
point, and are related to work done in BTl [38l l5l l26l [Tl . The reader should be warned that 
many of our proofs are quite technical in nature; we have tried to explain in advance the 
strategies behind the proofs, which are more conceptual. 

At the end of this paper we suggest several interesting questions related to the hypoel- 
lipticity of differential operators that are natural to consider in the context of distribution 
theory. It should also be noted that there are a number of results on local solvability of 
differential equations ll35l l24l that could be reformulated in this context. We expect that 
this work will provide the foundation for many subsequent investigations. 

2. Setting 

We begin by describing the basic elements of analysis on a post-critically finite self- 
similar set X, as laid out in the monograph of Kigami lTT~3l ; in this section all unreferenced 
results may be found in [13], which also includes proofs and references to the original 
literature. The reader who prefers to have a concrete example of a p.c.f. set in mind may 
choose to think of X as the Sierpinski Gasket, in which case an more elementary exposition 
of the material that follows may be found in |36|. 

P.C.F. Fractals. Let X be a self-similar subset of R d (or more generally a compact metric 
space). By this we mean that there are contractive similarities [Fj]^ of R d , and X is the 
unique compact set satisfying X = U^ =l Fj(X). Then X has a natural cell structure in which 
we associate to a word w - w\W2 ■ ■ ■ w m of length m the map F w = F Wl o • • ■ o F Wm , and call 
F W (X) an m-cell. If w is an infinite word then we let [w] m be its length m truncation and 
note that F W (X) = P| m F[ W ]JX) is a point in X. 

We say Fj(x) is a critical value of X = U^ =1 Fi(X) if there is y e X and k + j such 
that Fj(x) = Fk(y). An infinite word w is critical if F W (X) is a critical value, and w is 
post-critical if there is j e {1, . . . ,N] such that jw is critical. We always assume that the 
set of post-critical words is finite, in which case the fractal is said to be post-critically 
finite (p.c.f.). The boundary of X is then defined to be the finite set Vo consisting of all 
points Fa,(X) for which w is post-critical; this set is assumed to contain at least two points. 
We also let V„, = U w F w (Vb)> where the union is over all words of length m. Points in 
V* = U m >oV m that are not in Vo are called junction points, and a key property of p.c.f. 
fractals is that cells intersect only at junction points. 

We fix a probability measure p on X that is self-similar in the sense that there are 
pi, . . .,Hn such that the cell corresponding to w = W\ . . . w m has measure p(F w (X)) = 
n'jli fiwj- The usual Bernoulli measure in which each pj = jj is one example. 
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Dirichlet Form. Our analysis on X will be constructed from a self-similar Dirichlet form. 
A closed quadratic form E on L 2 (p) is called Dirichlet if it has the (Markov) property 
that if u e dom(£') then so is u — uxo <u <\ + Xu>\ and E(u,u) < E(u,u), where xa is 
the characteristic function of A. Self-similarity of E means that there are renormalization 
factors r\ , . . . , r^ such that 



It follows immediately that E(u, v) can also be expressed as the sum over m-words of 
r^E(u o F w , v o F w ) where r w = r W] ■ ■ ■ r Wm . In order to use results from [27 1 we assume 
that < rj < 1 for all j, in which case E is regular, meaning that that C(X) f] dom(£') 
is dense both in dom(£') with E'-norm and in the space of continuous functions C(X) with 
supremum norm. 

It is far from obvious that interesting fractals should support such Dirichlet forms, but 
in fact the conditions described so far are satisfied by many p.c.f. self-similar sets that have 
sufficient symmetry. In particular, if X is a nested fractal in the sense of Lindstr0m ifTTl 
then a Dirichlet form of the above type may be constructed using a diffusion or a harmonic 
structure 021 Eg). Some other approaches may be found in l20l[l4l [T8ll22ll2Tl[T0ll23l. 

Harmonic Functions. Given a function on Vq (usually thought of as an assignment of 
boundary values) there is a unique continuous function on X that has these boundary values 
and minimizes the energy. Such functions are called harmonic, and form a finite dimen- 
sional space containing the constants. It is easy to see that there are harmonic extension 
matrices Aj, j = 1, . . .,N with the property that if h is harmonic then Aj maps the values 
of h on Vq to its values on Fj(Vo). The largest eigenvalue of each Aj is 1, corresponding 
to the constant functions; it is useful to know that the second eigenvalue is rj, and that 
all other eigenvalues (which may be complex) have strictly smaller absolute value ( |[T3ll . 
Appendix A). 

The Laplacian and Normal Derivatives. Using the energy and measure we produce a 
weak Laplacian by defining / = Au if E(u, v) — — J fv d/i for all v £ dom(£') that vanish 
on Vo. Our assumptions so far are sufficient to conclude that -A is a non-negative self- 
adjoint operator on L 2 (fi) with compact resolvent (see Theorem 2.4.2 of fl~3)). We denote 
its eigenvalues by Aj and the corresponding eigenvectors by When Au e C(X) we 
write u e dom(A) and think of these as the (continuously) differentiable functions on X. 
Inductively define dom(A 4 ) for each k and then dom(A°°) = dom(A*). We say / is 
smooth if / E dom(A°°). Harmonic functions have zero Laplacian. 

By introducing a normal derivative d„ at boundary points the defining equation for the 
Laplacian can be extended to functions that do not vanish on Vo. As a result we have 
the Gauss-Green formula E(u, v) = - f(Au)vdn + 2xev v(x)d„u(x) when v £ dom(£), as 
in Theorem 3.7.8 of [13|. This formula may be localized to a cell F W (X), in which case 
d™u(q) - lim„, E(u, v m ) at the boundary point q = f] m F KI ^(X), where v m is the harmonic 
function on F wr m(X) with all boundary values equal to other than v m (q) = 1. The super- 
script w in d^u(q) indicates which cell the normal derivative is taken with respect to, as 
there is one for each cell that intersects at q. In general the normal derivatives exist once 
Am exists as a measure. If u £ dom(A) then the normal derivatives at a point sum to zero. 
Conversely, if u is defined piecewise by giving functions Uj £ dom(A) each supported on 
one of the cells that share the boundary point x, then u E dom(A) if and only if all Ui(x) are 




(2.1) 
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equal, all Aiij(x) are equal, and the normal derivatives of the Uj at x sum to zero. We call 
these constraints the matching conditions for the Laplacian. 

Resistance Metric. In addition to the Laplacian and other derviatives, the Dirichlet form 
also provides us with a metric intrinsic to the fractal. We define the resistance metric R by 

R(x,y) = tmn{E(u) : u e dom(£'), u(x) = 0, u(y) — 1}. 

In Sections 2.3 and 3.3 of [13] it is proven that under our assumptions this minimum exists 
and defines a metric, and that the /^-topology coincides with the topology induced from the 
embedding of X into R , Of particular importance for us is the fact that continuity may be 
treated using the resistance metric, for which purpose the following Holder- i estimate is 
very useful: 

\u(x) - u(y)\ 2 < E(u)R(x,y) < ||ii|| 2 ||Aii|b*(JC,;y). (2.2) 

If u € dom(A) vanishes on Vo then we obtain the first inequality trivially from the definition 
and the second by applying the Cauchy-Schwartz inequality to E(u) - - J uAud/j. For 
general u e dom A we can simply subtract the harmonic function with the same boundary 
values and apply the same estimate. In particular, this shows that the Lr domain of the 
Laplacian embeds in the continuous functions. 

Fractafolds. Since the results in this paper are primarily local in nature, we will be able to 
work on a connected fractafold based on X with a restricted cellular construction, which we 
denote by SI. Some results on fractafolds and their spectra may be found in [32). As with a 
manifold based on Euclidean space, a fractafold based on X is just a connected Hausdorff 
space in which each point has a neighborhood homeomorphic to a neighborhood of a point 
in X. One way to construct a fractafold is by suitably gluing together copies of X, for 
example by identifying appropriate boundary points. This leads us to the idea of a cellular 
construction, which is the analogue of a triangulation of a manifold. A restricted cellular 
construction consists of a finite or countably infinite collection of copies Xj of X, together 
with an admissible identification of their boundary points. Admissibility expresses the 
requirement that the result of the gluing be a fractafold; more precisely, it means that if 
\x\, . . . ,xj] are identified then there is a junction point x e X and a neighborhood U of x 
such that each of the components U\, . . . , Uj of U\ {x} is homeomorphic to a neighborhood 
of the corresponding point xj in Xj, We call any such point x a gluing point, and make the 
obvious definition that a neighborhood of x is a union of neighborhoods of x in each of the 
cells Xj that meet at x in the manner previously described. 

It should be noted that the above is not the most general kind of cellular construction 
(hence the term restricted in the definition), because some fractals have non-boundary 
points (called terminal points) at which cells may be glued (see [32], Section 2). Dealing 
with such points introduces certain technicalities that, while not insurmountable, cause 
complications in defining the Green's operator (see below) that we will need for proving 
Theorem |5.9l It is worth noting that if X has some topological rigidity then all fractafolds 
have restricted cellular structure. This is true, for example, for fractafolds based on the 
Sierpinski Gasket ( 021 Theorem 2.1). 

Thus far our fractafold has only topological structure; however if SI has a restricted 
cellular construction then a smooth structure may be introduced in the same manner as it 
was on X itself, specifically by defining a Dirichlet energy and a measure and thus a weak 
Laplacian. We can take the energy on SI to be the sum of the energies on the cells of 
the cellular construction, and the measure (which is not necessarily finite, but is finite on 
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compacta) to be the sum of the measures on the cells: 

E(u, v) = E Xj (u\ x . , v\ x ) = a J E M Xj ° l p v \x, ° h) 
j j 

J 7 

where ij : X —> Xj is the map from the cellular construction. In the same way that the angle 
sum at a vertex of a triangulation of a manifold determines the curvature at the vertex, the 
choice of the weights aj and bj amount to a choice of metric on Q, with cij — bj = 1 for 
all j being the flat case (see (32), Section 6). As all of the computations made later in the 
paper may be made on one cell at a time, we will henceforth suppress the weights cij and 
b, 

Well-known results about p.c.f. fractals imply the existence of a Green's function (for 
which there is an explicit formula) on finite unions of cells in a fractafold with cellular 
construction. 

Lemma 2.1. Let K — u|Xy be a connected finite union of cells in Q and such that K + 
Q. Then there is a Green 's operator Gk with the property that if v is a Radon measure 
on K (i.e. a Borel measure that is finite on compacta, outer regular on Borel sets and 
inner regular on open sets), then Gkv is continuous, —AGkv — v on the interior of K, 
and GKV\ gK — 0. The same conclusion holds in the case K — Q under the additional 
assumption J dv — 0. 

Remark 2.2. It is clear that dK is a subset of the boundary points of the cells Xj, specifi- 
cally consisting of those gluing points at which not all glued cells are included in K. 

Proof. We recall from Sections 3.6 and 3.7 of 0~3| that our assumptions on X imply there 
is a Green's operator G on X with continuous kernel g(x,y), such that -AGv - v and 
Gv| ex = for all Radon measures v. There is an explicit formula giving g(x,y) as a series. 

If Gj is the Green's operator for the cell Xj it is easy to verify that -A £ GjV = v, 
except at the gluing points where the Laplacian can differ from v by Dirac masses, the 
size of which may be computed explicitly by summing the normal derivatives of the GjV 
at the points that are glued. However it is also apparent that by assigning values at each 
of the gluing points and extending harmonically on the cells we obtain a continuous and 
piecewise harmonic function, the Laplacian of which is a sum of Dirac masses at the gluing 
points. 

Provided the boundary dK is non-empty (which is obvious if K + Q), a linear algebra 
argument (Lemma 3.5.1 in |[T3l ) shows that for any prescribed set of weights for Dirac 
masses of the Laplacian at interior gluing points, there is a unique piecewise harmonic 
function that vanishes on dK and has this Laplacian. Subtracting this piecewise harmonic 
function from 2 GjV gives the required Gkv. 

On the other hand, if K = O then Q is compact and the kernel of A is precisely the 
constant functions. We can therefore invert -A on the measures that annihilate constants, 
that is, those for which J dv — 0. This can be done explicitly in the same manner as in the 
previous case, except that the linear algebra step now shows the Laplacians of piecewise 
harmonics span the space of mean-zero linear combinations of Dirac masses at the gluing 
points. In this case the choice of piecewise harmonic function is unique up to adding a 
constant; our convention is to choose this constant so J G/j-v(x)t//j(x) = 0. □ 

Throughout this paper we will assume that Q. has no boundary. In some examples it is 
possible to deal with boundary points by passing to an appropriate cover, but relatively little 
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is known in terms of covering theory for general fractafolds. Elementary examples to keep 
in mind include non-compact cases like open subsets of X \ Vo or the infinite Sierpinski 
Gasket treated in 071 . and compact fractafolds like the double cover of the Sierpinski 
Gasket SG, which consists of two copies of SG with each boundary point from one copy 
identified with exactly one boundary point of the other (see [32] for more details). 

Smooth Cutoffs and Partitioning. As mentioned earlier, the structure theorems we shall 
prove for distributions rest heavily on results from 11271 . In what follows we assume that 
£2 is a fractafold with restricted cellular structure and is based on a fractal X with regular 
harmonic structure. 

Recall that x e X is a junction point if and only if there is a neighborhood U 3 X such 
that U \ \x] is disconnected into a finite number of components Uj. For a smooth function u 
the quantities A k u(q) and d J „A k u{q) exist for all k e N; the superscript / on d J „ indicates the 
normal derivative with respect to the cell Uj. For a fixed j, the two sequences A k u(q) and 
d J n A k u(q) make up the jet of u at q in Uj. The first result we need from ll27l is a Borel-type 
theorem on the existence of smooth functions with prescribed jets. 

Theorem 2.3 (J27j, Theorem 4.3 and Equation 4.8). Given values po,Pi, ■ ■ ■ ando-Q, <t\, . . . 
there is a smooth function f on Uj that vanishes in a neighborhood of all boundary points 
except q, where the jet is given by A k f(q) — pi and d' n A k f(q) — o~kfor all k. If we write Uj 
as Uj — F W (X) for a word w, and fix a number L of jet terms, then for any e > we may 
construct f so that for < k < L, we have the estimate 

L L-l 

||A*/IU < CftXr^nJ r'J w \ P i\ + Yj # V>/l) + c (2-3) 

1=0 1=0 

where C(k) depends only on k and the harmonic structure on X. 

Remark 2.4. Of course, it follows immediately that we can construct a smooth function 
with prescribed jets at each of the boundary points of a cell K and an estimate like ( 12.31 l. 
just by applying the theorem separately to each of the boundary points and summing the 
result. 

Corollary 2.5. // K is a cell in Q and U is an open neighborhood of K, then there is a 
smooth function f such that f — 1 on K, f — outside U, and ||/||oo < C, where C is a 
constant that does not depend on K or U. 

Proof. Let {qj} be the boundary points of K and at each qj take cells Vjji c U such that 
U/fc Vjji u K contains a neighborhood of qj. By making all of these cells sufficiently small 
and removing any inside K we may further assume that the Vjjc intersect K only at qj, 
intersect each other only at qj and do not intersect Vfj? for any / + j. 

On each V}^ construct the smooth function fj^ guaranteed by Theorem l2 . 3 1 with = 1 
at qj and all other jet terms at q equal to zero, and taking e = 1 . Then the piecewise 
function 





1 


for x e K 


/(*) = • 


fj* 


for x 6 Vj 







otherwise 



is equal 1 on K and off U by construction. It is also smooth, simply because the pieces 
are smooth and the matching conditions for A' apply at each of the boundary points of the 
Vj,k for all I. The bound H/IU < C independent of K and U now follows from j2.3i because 
the scale-dependent terms are all raised to the power zero, so are constant. □ 
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A more difficult task than that in Corollary 12. 5l is to construct a positive bump function 
that is equal to 1 on A" and to zero outside the neighborhood U of K. A result of this 
type was proven in 11271 under certain assumptions on the diffusion Y, corresponding to the 
Laplacian. A sufficient assumption is that the heat kernel p,(x, y) (i.e. the transition density 
of of Y,) satisfies an estimate of the form 

where a, fl, y\ and j2 are constants. The estimate [2.4| is known to be valid in great gener- 
ality on p.c.f. fractals (Corollary 1.2 of [9|). 

Theorem 2.6 ( 1271 Corollary 2.9). Under the assumption ( 12.4| >, for a cell K and an open 
neighborhood U D K, there is a smooth function f such that f — 1 on K, f — outside U, 
and f{x) > for all x. 

The final theorem from [27] that we will use extensively is concerned with the smooth 
partitioning of a smooth function. 

Theorem 2.7 ( 11271 , Theorem 5.1). Let K c X be compact and fix \J U a "J K an open 
cover. If f e dom(A°°) then there is a decomposition f = Y/ifj m which each fj is in 
dom(A°°) and has support in some U a .. 

Remark 2.8. Compactness of K is used only to obtain finiteness of the decomposition, 
and may be omitted for finite covers. An analogous countable (and locally finite) decom- 
position is then valid in the cr-compact case; in particular it is valid on Q, because of the 
existence of a cellular structure. 

Remark 2.9. The proof uses a result on the existence of smooth functions with prescribed 
jet at a point (Theorem 4.3 of [27]) to smoothly join cutoffs to a piece of the original func- 
tion as in the proof of Theorem 12. 5 1 This is very different from the Euclidean case where 
one simply multiplies the smooth function by a smooth bump. In particular, the construc- 
tion of the cutoff depends explicitly on the growth rate of the jet of / at the boundary points 
under consideration, so for a collection of sets indexed by j, the mapping / i-» [fj] to a 
sequence of smooth functions supported on these sets is nonlinear. 

Although the non-linearity will make some later proofs more complicated, this method 
does provide good estimates. From ( 12.31 ) and standard arguments for controlling the normal 
derivative d„A k f at a point by the norms HA^U, j = . . . , k + 1, over a neighborhood of 
the point (like those in Section|6]below) we find that / h-> fj can be arranged to satisfy 

k 

IIAV-IU <C^J|A'/IU (2.5) 

1=0 

where C is a constant depending only on k and K. 

3. Test Functions 

We define test functions on Q. in the usual way, and provide notation for the space of 
smooth functions on Q. topologized by uniform convergence on compacta. 

Deflntion 3.1. The space of test functions £)(Q) consists of all <\> e dom(A°°) such that 
Sppt(0) is compact. We endow it with the topology in which 0,- — > iff there is a compact 
set K c Q containing the supports of all the </>,, and A*0, — > A k (p uniformly on K for each 
teN. There is a corresponding family of seminorms defined by 

|0L = sup{|AVW| : x € £l,k < m) (3.1) 
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though it should be noted that the topology on TXQ.) is not the usual metric topology 
produced by this family. For a discussion of the topology on £)(£!) and its relation to these 
seminorms, see Chapter 6 of [28 1. 

Deflntion 3.2. fi(Q) = dom(A°°) with the topology <p-, — » cf> iff for every compact K c Q 
we have A k (f>i — » A k <p uniformly on K for each k eN. There is a corresponding family of 
seminorms defined by 

\<P\ mX = sup{|AVM| :xeK,k< m}. (3.2) 

The following result is immediate from Theorem l2.7l and ( 12.5b . It will be used frequently 
in the results proved below. 

Lemma 3.3. If <p e 2)(Qi U Q.2), then <p = <p\ + $2 for some <pj e D(Q.j). For each m there 
is C = C(M, a u Q 2 ) so \<Pj\ m < C|0L 7 = 1,2. 

One consequence is that D(Q) is dense in £(£2), because we may fix an increasing 
compact exhaustion UjKj = Q of our domain and for arbitrary <p e S(Q.) write <p — <pj + 
where d> t is supported in K,+\ and (p ,■ is supported in /f c , so that = ; . The functions 

J J J j \ K j ■> \ Kj 

<pj are in D(Q) and it is clear that A k <f>j — > uniformly on compacta, hence — > in 
Another density result that follows from Lemma [331 is as follows. 

Theorem 3.4. £)(Q) is dense in C C (Q), the space of continuous functions with compact 
support, with supremum norm. 

Proof. The dual of C C (Q) is the space of Radon measures, so by the Hahn-Banach Theo- 
rem, it suffices to show that if such a measure v satisfies 



dv = 0, for all E £>(Q), (3.3) 



then veO. 

Let K be a cell and {t/, } a sequence of open sets containing K so that v(t/, \ K) — > 0. 
Using Corollary 12.51 we see that for each j we can take 0, e 2D(f2) with 0, = 1 on /f, the 
bound ||0,||oo < C for all i, and Sppt(0,) c (!/;)• Then for v satisfying (13 - 31 > we compute 



y(K) 



s v(^ \ X) < Cv(f/j \ X) 0. 



As v vanishes on all cells it is the zero measure, and the result follows. □ 

Since Q. is locally compact and Hausdorff, it is a standard result that C f (Q) is supremum- 
norm dense in Co(£2), where the latter consists of those continuous functions / for which 
the set {x : \f(x)\ > e] is compact for all e > 0. Hence D(Q) is also dense in Cq(Q). 

In the special case where Q. is compact we may also characterize D(Q) = S(Q.) by 
the decay of the Fourier coefficients obtained when <p is written with respect to a basis 
of Laplacian eigenfunctions. This provides an alternate proof of the density of £)(Q) in 
C C (Q), which of course coincides with C(£2) in this case. 

Theorem 3.5. If £2 is compact then £)(Q) = is the space of smooth functions with 
Fourier coefficients that have faster than polynomial decay, and hence is dense in C(f2). 

Proof. Clearly (f> e D{Q) is in L 2 , so can be written <p = Yj^o a i4 1 ^ where fa is the Laplacian 
eigenfunction with eigenvalue -A/. It follows that (-A) k <p = 2/ diAtyi with convergence 
in L 2 . Since A k <p is in C(Q.) c L 2 for all k we see that the sequence a, must decay faster 
than any polynomial in the i,. Conversely any such sequence converges to a function 
for which every power of the Laplacian is in L 2 , whereupon the function is smooth by 
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iteration of (12.2b . In addition, any u e C(Q) can be explicitly approximated by functions 
from D(O) by taking successive truncations of the Fourier series u = £°! a^,. To see this 
gives convergence in D(£2) write (—A)* a i t / / i - lLT=j a i^\^i ar, d note this converges to 
zero in L 2 and therefore almost everywhere. Now from d2.2| ) 

CO ^ CO CO 

su P |2 «/4<A,| < c|J] fl^V|J|Z fl 4 + Vi|| 2 

and both terms on the right converge to zero. □ 

In ll33ll there is a definition of Sobolev spaces on p.c.f. fractals of the type studied here. 
These spaces may be defined by applying the Bessel potential (/ - A)~ s (for the Dirichlet 
or Neumann Laplacian) or the Riesz potential (-A) s (for the Dirichlet Laplacian) to LP 
functions on the fractal, and adding on an appropriate space of harmonic functions. In 
particular, the space of L? functions with A k u e L? for < k < m may be identified with 
a particular L? Sobolev space ([33] Theorem 3.7). Writing W s ' 2 for the L? Sobolev space 
arising from (/ - A)~ s , we have in consequence of the preceding: 

Corollary 3.6. If Q is compact then £>(Q) = n s>0 W s ' 2 . 



4. Distributions 

Deflntion 4.1. The space of distributions on Q. is the dual space £)'(Q) of D(£2) with the 
weak-star topology, so T, — > T if and only if — > T(p for all (p e D(Q). 

As usual, the most familiar examples of distributions are the Radon measures. If dv is 
such a measure then we define T v by T v <p = J cf>dv. Theorem l3.4l ensures that the mapping 
y i — > T Y is injective, so we may identify v and T v . One way to obtain further examples 
is to take the adjoint of the Laplacian on distributions, which clearly produces another 
distribution. 

Deflntion 4.2. If T € £)'(Q) we define AT e £>'(Q) by (AT)<f> = T(A<f>) for all <p e 0(Q). 

It is clear that powers of the Laplacian applied to the Radon measures provide a rich col- 
lection of examples of distributions. Later we prove that all distributions arise in essentially 
this way (Theorem |5.9l ), but we first need to establish some more elementary properties. 

Theorem 4.3. A linear functional T on £)(Q) is a distribution if and only if for each 
compact K c Q there are m and M such that 

\T(f>\ < M\(f>\ m (4.1) 

Proof. It is clear that the existence of such an estimate ensures continuity of T. To prove 
the converse we assume no such estimate exists, so there is K compact and a sequence <f>j 
such that \T<f>j\ > j\<f>j\j. Then the support of <f>j = 4>jlT4>j is in K for all j and 

\\A k d>j\\ < —— < - once / > k. 
II ^IL \ T( f>j\ j J 

Therefore <pj — > in D(Q) but has T<pj = 1 for all j, contradicting the continuity of T. □ 

In the special case that Q is compact we saw in Theorem 13.51 that D{Q) consists of 
smooth functions having Fourier coefficients that decay faster than polynomially. This 
allows us to identify its dual with coefficient sequences having at most polynomial growth. 
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Lemma 4.4. IfD. is compact and T e D'(Q) then the sequence Ttffi has at most polynomial 
growth. Conversely, any sequence {£>,} of polynomial growth defines a distribution via 
<f> = ZiOiif/i YjiUibi. 

Proof. We saw in Theorem [33] that Ylt = n a 0l'i converges to in D{Q) if and only if {a,} 
has faster than polynomial decay in A/. It follows that for any T e D'{Q), a{T^ii = 

T YJi = Q Qfl/fi ~ * T<p, from which the sequence Txji; has at most polynomial growth. 

Conversely suppose that {/?,} has polynomial growth, < CA m , and consider the map 
{fl/} Tii a ibi- This is a well defined linear map T on sequences {a,} with faster than 
polynomial decay, hence on £)(Q), with the estimate 

17-01 < J] \aM < Cj] W = C||A'>|| 2 < C||A'"0|| sup < C|0L 

i i 

which shows that T is a distribution. □ 

In particular, if we identify £)(Q) as a subset of £)'(Q) by letting 0' e £>(Q) act on £)(Q) 
via h-> (0, 0'), where (, ) is the L 2 inner product, then this implies that the test functions 
are dense in the distributions when Q. is compact. To see this, define Tj by 



|7>; ifi<j 
[0 if»>/ 

We see that (T - Tj)(p = a/b, —* for any (f> E £)(Q), so 7^ — > 7*. Since Tj is the inner 
product with the function idf>u it is in This is true more generally. 

Theorem 4.5. £)(Q) is dense in 

Proof. Let T e £)'(Q). Take an increasing exhaustion UKj of f2 by compact sets Kj with 
the property that Kj is contained in the interior of Kj+\, and each Kj is a finite union of 
cells. For each j, the action of T on DiKj) identifies it as an element of D'(Kj) so by the 
preceding there is a sequence {Tj^T-n c 2)(Kj) for which 7y^ — » T in D'(Kj), and hence 
in £)'(£/) for all / < j. 

Now consider the diagonal sequence Tj j. For any test function (f> there is some jo such 
that Kj contains the support of 0, and hence Tjj(/> is defined for j > jo and converges to 
T(f>. So Tjj — » r in Of course, 7% ( y only corresponds to a test function ; on Kj, not 

to an element of £)(Q). To remedy this, note that for the test function <pj corresponding to 
Tjj on Kj we may apply Theorem l2.3l to each of the (finite number of) boundary points of 
Kj and thereby continue <pj smoothly to a function (p'j on Q. with support in Kj+i. Denote 
by Tj the distribution corresponding to this new test function 4>'j- Since 4>j and <p'. coincide 
on Kj we see that T )</> = T'.<j> for all <p having support in Kj. It follows that T. converges to 
T in 2)'(Q), and since each T. corresponds to a test function, the proof is complete. □ 

Deflntion 4.6. If Q.\ cQ is open, we say the distribution T vanishes on Qi if Tip — for 
all <p supported on Cli, This is written T\ n =0, 

To make a meaningful definition of the support of a distribution we again need the 
partitioning property. By Lemma [331 we know that any <j> € D(Q.\ U Q2) can be written as 
4> = 0i + <f> 2 for 4>j € £>(Q ; ). If both Tfr = and T0 2 = it follows that Ttf> = 0. We 
record this as a lemma, and note that it establishes the existence of a maximal open set on 
which T vanishes. 

Lemma 4.7. 7/tL =0andTL =0thenTL n = 0. 



DISTRIBUTION THEORY ON P.C.F. FRACTALS 



11 



Defintion 4.8. The support of T is the complement of the maximal open set on which T 
vanishes, and is denoted Sppt(T). In the special case where Sppt(T) is compact we call T 
a compactly supported distribution. 

Theorem 4.9. The space of compactly supported distributions is ( naturally isomorphic to ) 
the dual £'(£2) of the smooth functions on Q. 

Proof. The inclusion £)(Q) c £(Q) defines a natural map from £'(Q) to £)'(Q). We have 
seen (after Lemma [33] above) that 2D(Q) is dense in £(Q), from which it follows that the 
kernel of this map is trivial. For convenience we identify £'(£2) with its isomorphic image 
in £)'(Q), so we need only verify it is the space of distributions with compact support. 

Fix an increasing sequence of compacta Kj with Q = UjKj. If T € £'(Q) fails to be 
compactly supported then for each j there is <f>j supported in Q. \ Kj such that T<f>j + 0, and 
by renormalizing we may assume T(pj = 1 for all j. However for any compact K there is 
j such that K c Kj and thus 0/ = on K once I > j. This implies <pj — > in £(£2) and 
T<pj = 1 for all j, a contradiction. 

Conversely, let T e D'(Q) be supported on the compact K c Q. Fix a strictly larger 
compact K\ (so that K is contained in the interior of K\) and an open neighborhood Oi of 
K\. By Remark l2~8l the conclusion of Theorem l2.7l is valid for the cover by Q. \ Ki and Qi, 
even though £2 is noncompact. In particular we can fix a decomposition mapping in which 
/ e £(£2) is written as / = f\ + f%, with fx supported on Q. \ K\ and therefore f\ | = /| . 
Now let T\ on £(Q) be given by T\f = Tf\. This is well defined, because if /, g € £(Q) 
and /] = gi, then f — g = in a neighborhood of K and the support condition ensures 
Tf = Tg. It is also linear, even though the mapping / \-> fx is nonlinear (see Remark 
12.9b . because (/ + g)\ - f\ + gi on K\, which contains a neighborhood of K. Lastly, T\ is 
continuous, as may be seen from the fact that a sequence {<f>j} c £(Q) such that A k (f>j — > 
uniformly on compacta will have A k {(pj)i — > on K\ D K, or from ( 12.5b . We conclude that 
every compactly supported distribution is in £'(f2). □ 



5. Structure Theory 

Defintion 5.1. A distribution T has finite order m if for each compact K there is M = M(K) 
such that \T(f>\ < M\<f>\ m for all cf> € £>(#). 

The following theorem indicates the importance of the finite order distributions. 

Theorem 5.2. Compactly supported distributions have finite order. 

Proof. Let T be a distribution with compact support K and let K\ be a compact set such 
that K c mt{K\). By Lemma l^Ol we may decompose any 6 £)(Q) as = <pi + 02 where 
0i is supported on Xi, 02 is supported in A^', and \<f>j\ m < C m \<p\ m for j = 1,2. Clearly 
T0 = T0i , but there are m and M so that ( 14. U holds on /Tj , from which we conclude that 

|T0| = |7>i| < M|0iU < C m M|0| m . □ 

It is easy to see that the Radon measures on £2 are examples of distributions of finite 
order. In fact they have order zero, because the action of v on £)(Q) via v0 = J if>v trivially 
satisfies the bound |v0| < H0IU = |0lo- The converse is also true. 

Theorem 5.3. If T is a distribution of order zero then there is a Radon measure v such 
that T<p = f <f>dvfor all £ £)(Q). 
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Proof. Let K be compact. Since T has order zero there is M — M(K) so that \T<p\ < M\\<p\\oo 
whenever e D{Q) has support in K. Theorem B. 41 shows that these functions are dense 
in C(K), so we may extend T to a bounded linear operator on C(K). Such operators are 
represented by Radon measures, so there is vk with T<p — J fidvx for all test functions (p 
with support in K. Now let |J Kj be a compact exhaustion of Q and consider the measures 
Vk ■ These converge weak-star as elements of the dual of C e (Q) to a Radon measure v, but 
by construction J 4>dvK t — > T<p on D(£2) and the result follows. □ 

Remark 5.4. As written, the preceding proof relies on Theorem 13.41 and hence on the 
Hahn-Banch theorem. Since each of the Kj is compact, a constructive proof can be ob- 
tained by instead using Theorem l3~5l 

A well known application of the preceding is obtaining a characterization of the distri- 
butions that have positive values on positive test functions. To prove this we need Corollary 
12.61 and therefore must make the corresponding assumption ( 12.41 ) on the behavior of the 
heat kernel corresponding to the Laplacian. 

Deflntion 5.5. T is a positive distribution if T(f> > whenever <p > is a positive test 
function. 

Theorem 5.6. Positive distributions have order zero. If the Laplacian on X is such that 
(12.4b holds and if T is a positive distribution, then there is a positive measure v such that 
I f - ffdv. 

Proof. Let K be compact. Using Theorem 12.51 there is if/g e D(Q) such that iff = 1 on A". 
If (p e C°° with support in K then the functions H^U^x ± 4> are both positive, whence 

-MUTif/K <T4>< U\Uiff K . 

We conclude that T has order zero, so by Corollary 15.31 it is represented by integration 
against a measure v. If there is a cell K for which v{K) < then we can take Uj to be a 
neighborhood of K for which v(Uj \K) < l/j and let fj be as in Theorem l2.6l It follows 
that 

Tfj = Jfdv<v(K) + j\\f\L 

and for a sufficiently large j this is negative, in contradiction to the positivity of T. We 
conclude that v(K) > for all cells K, and therefore that v is a positive measure. □ 

We noted at the beginning of Section|4]that the adjoint of the Laplacian maps £)'(0) to 
itself. In particular, if v is a Radon measure, hence a distribution of order zero, then for 
each compact K there is M{K) such that 

|(A'"v)0| = |v(A'»| < M(K)\A m <f>\ < M(K)\(p\ m 

so A m v is a distribution of order m. This result has a converse, which we prove using a 
modification of the Green's function introduced in Lemma lXTI The basic idea is to produce 
a Green's operator that inverts the Laplacian on test functions, so that the adjoint of this 
operator lowers the order of a finite-order distribution. Iterating to produce a distribution 
of zero order then produces a measure by Theorem |5.3l 

Lemma 5.7. Let K be a connected finite union of cells in Q. Then A : D(K) — > D(K) and 
its image consists of all test functions that are orthogonal to the harmonic functions on K. 
Moreover there is a linear operator Gk '■ IXK) — > D(K) such that —AGk(A(P) — A(f> for 
all 4> e D(K). 
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Proof. For <p e D(K), the matching conditions for the Laplacian ensure that both (p and d„(p 
vanish on dK because <p is identically zero outside K. If/ e dom(A), then the Gauss-Green 
formula reduces to 



f (Acf>)fdfi = f 0(A/)4u 



because there are no non-zero boundary terms. We conclude that / is orthogonal to the 
image A(D(K)) if and only if A/ is orthogonal to D(K). As the latter is dense in C(K) 
(Lemma [3. A\ the first result is proven. 

Let h\ . . . , hi(iq be an L 2 -orthonormal basis for the finite dimensional space of harmonic 
functions on K. As D(K) is dense in C(K) in supremum norm and K has finite measure, 
D(K) is also dense in both C(K) and L 2 (K) in L 2 norm. It follows that there are (f>i, . . . , <Phk) 
in D(K) such that {<pi,hj) = where (, ) is the L 2 inner product and <5,y is Kronecker's 
delta. Given <f> e D(K) we let 

m 

4> = Y J (<f>,h i )cf> i (5.1) 



and define 

G K 4>(x) = 



G K (<p-<f>)(x) ifxeK 
if jc ^ A" 

where Ga: is the Green's operator defined in Lemma |2TI It is then clear that for if/ e D(K), 

- AGki/s = -AG K (i/f -$)=\{/-i, (5.2) 

except perhaps at points of dK, where we must first verify that the matching conditions for 
the Laplacian hold. Since Gk^ vanishes outside K, the matching conditions require that 
dnGK^iq) = whenever q e dK. One way to verify this is from the Gauss-Green formula 
for a harmonic function h, which yields 

Q = ty-4r,h)= f (-AG K iff)h = ~Yj (-d n G K Hq))Kq) 

^ K qedK 

from which we see that it suffices to know the solvability of the Dirichlet problem on K, 
that is, for every assignment of boundary values on dK there is a harmonic function h 
with those boundary values. This latter is true because of Lemma |2~T1 for example it may 
be proven by taking a function that is piecewise harmonic on cells and has the desired 
boundary data and subtracting the result of applying Gk to its Laplacian (which is simply 
a sum of Dirac masses at the interior gluing points). We conclude that Gk^ 6 35 (K) and 
that ( 15.2b holds everywhere. 

Finally, if if/ — Acp for some <f> e T>(K), then if/ is orthogonal to the harmonics, so if/ = 
and -AGk^ — iff as desired. □ 

The adjoint of Gk is defined on distributions by {GkT)4> = T{Gk4>)- This operator is 
really defined on the dual of D(K), which is a larger space, but we will not make use of 
this fact. 

Theorem 5.8. If T is a distribution of order m > 1 then GkJ is a distribution of order 
m — 1, and if T is a distribution of order zero then GkT is integration with respect to a 
continuous function on K. 

Proof. Let T be a distribution of order m > 1, so that GkT<P is bounded by 
\G K <p\ m = supfllA^G^IL : k < m). 
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When k > 1 we have A k G K {<!>) = -A k ~ l ((f> - 4>), and when k = we see that ||Gjc(0)||so < 
C\\(f> - 4>\\<x, because the operator Gk in Lemma |2~T1 is clearly bounded on L°° . Hence 
\Gk<P\ih - C\<p - 4>\ m -\ < C(m - \,K)\<f>\ m -i, where the latter inequality is from ( 15. Q with 
a constant C(m -\,K) that may depend on the set of functions <pj. Thus GkT has order 
m — 1. 

If T has order zero then by Theorem l5.3l it is represented by integration against a Radon 
measure v. Provided K + Q. we can apply Lemma I27TI directly to see v = Af for some / 
that is continuous on K and vanishes on dK, so can be extended continuously to be zero 
outside K. This ensures there are no boundary terms when we compute with the Gauss- 
Green formula: 

G K T<p = TG K (f>= f G K <pdv= f (G K <f>)(Af)dfi 
Jk Jk 

= f (-AG K cf>)fd»= [ to-fofdji 
Jk Jk 

= I (0-^<<Mi>&)/4" 



!=1 



{(f>i,f)hi)dfi 

i=i 

and the bracketed term in the last line is continuous because / is continuous and all of the 
hi are harmonic. 

The argument is slightly different if O = K. We instead set t = J dv/(J d/j) so J d(v - 
fyu) = 0, at which point Lemma lZTI applies to show v-tfi- Af for a continuous /, and we 
can compute as before 

GkTc/) — I Gn<pdv 
Jk 

= t f G K 4>d[i+ f G K (f>(Af)dfx 
Jk Jk 

(<f> - cf>) f dn 

c 

= [ 4>{f- f fdix)dn 

where we used that f Gk4> dfx = (from the proof of Lemma 12. U and that the harmonic 
functions are constants in this case. □ 

We now have all the necessary tools to prove the main structure theorem for distribu- 
tions. 

Theorem 5.9. Any distribution T may be written as a locally finite sum of the form T — 
2 A m 'Vj or T — J] & mj+l fp where the Vj are Radon measures and the fj are continuous 
functions with compact support. 

Proof. Suppose first that O is non-compact and take K\ , K2, ... a sequence of subsets such 
that each Kj is a connected finite union of cells, Kj is contained in the interior of Kj+\, and 
*JjKj = Q. Such a sequence exists because Q. has a restricted cellular construction. It will 
be convenient to also set Kq = 0. For each j let Gj = Gk, be the operator from Lemma lBTTl 
The key point of the proof is that for any distribution S , we have (-A) m G'JS — S as 
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elements of D'(Kj) (though not as elements of D'(Q)). This may be verified by direct 
computation. For all <p e D(Kj), 

- AGjScf> = (-GjS)(A(f>) = S(-Gj(A<f>)) = -S<p (5.3) 

where the final step uses that -AGj(A<p) = A</> from Lemma l5"771 so G/(A0) + (f> is harmonic 
on K and vanishes on dK, hence is identically zero. 

Fix T e £)'(Q). Inductively suppose that for i = 0, ... ,j — 1 we have found m, and a 
measure v, supported on K, such that T - £jT A" 1; v, vanishes on D(Kj-i). The base case 

j = is trivial because A"o = 0. Now T - 2q 1 A m, v; is in £>'(£!), hence its restriction to 
£)(Kj) is in D'{Kj). We call this restriction Tj. As is compact, Tj has finite order »iy. It 
satisfies Tj = (-A) m 'G"j'T j by the argument already given, and by Theorem 15.81 there is a 
measure Vj supported on Kj such that Vj = (—Yfj&. J Tj. Therefore Tj = A'"'Vj in D'(Kj), 

which is equivalent to saying that T ~YJq A" 1 ' v, vanishes on D{Kj). 

It is immediate from the definition that £ ; A m 'vy is a locally finite sum. If we fix e 
£>(Q) then there is a j such that e D(Kj), whereupon (T - ^- =1 A'"'v,)0 = for all I > j. 
This proves that T = 2# A m 'vj. 

The proof that T = Yij & m ' +l fj is similar. Obviously we wish to use the latter part of 
Theorem 15. 8 1 to go from the measure vj to a continuous function. The only technicality 
is that the resulting fj is continuous on Kj rather than on all of Q. We fix this at each 
step of the induction as follows. Suppose we have determined Tj as the restriction of 
T - Tj\Zq A m,+1 fj and from Theorem l5.8l a function gj continuous on Kj such that A' n > +1 gj = 
Tj in D'{Kj). Let fj be a continuous extension of g ; to O obtained by requiring fj = 
on dKj + i and outside Kj+\, and letting it be piecewise harmonic on the cells of the cellular 
structure on Kj+\ \ Kj (here we use that Kj is in the interior of Kj + \). Clearly A mj+1 fj = Tj 
in D'(Kj), because we fj = gj on Kj, so T - % J i=0 A m, * +1 /i vanishes on !D(Kj) and we may 
complete the proof as before. 

In the case when Q is compact the proof is somewhat more elementary because we need 
only a single set K = Q, but there is a small technical difference due to the fact that the final 
equality of ( 15.3b is no longer true. Indeed, Gj(A(p) + <p is harmonic by the same reasoning as 
for the non-compact case, but now it is the possibly non-zero constant 4> = (ji(K))~ x J (f> d/j. 
The analogue of ( 15.3b is therefore 

AG K S<p = S<p-$Sl 

where 1 is the constant function. 

The distribution T has finite order to, and T(f> = Tl(f> + AGxT<p. Iterating, we have 
G K T(p = G K T10 + AG 2 K Tcj), and then 

L 

T<l> = A L G\T<l> + \Ya3 
1=1 

where a/ = & K T1. If L - to then Theorem 15.81 implies the first term is A"'v, where v is a 
measure, and if L — m + 1 this term is A m+l f, where / is a continuous function. In either 
case the second term is a constant multiple of the measure //, or equivalently the constant 
(hence continuous) function 1, so the proof is complete. □ 

6. Distributions supported at a point 

A distribution with support a point q is of finite order by Theorem 15.21 and simple 
modifications of the arguments in Theorem 15.91 show that it is a power of the Laplacian 
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applied to a measure with support in a neighborhood of q. The purpose of this section is 
to identify it more precisely as a finite sum of certain derivatives of the Dirac mass at q; in 
general these derivatives are not just powers of the Laplacian, but instead reflect the local 
structure of harmonic functions at q. 

Identification of a distribution T of order m supported at q is achieved by describing a 
finite number of distributions Tj,j- 1 , . . . , J with the following property: if <p e T)(Q) has 
Tj(p = for all j then for any e > there is a neighborhood U e of q and a decomposition 
<p = <p q + (<f> — <p q ) into test functions such that \<p q \,„ < e and <p - <p q vanishes on U € . The 
reason is that then T<p - T<p q because of the support condition, and \T<p q \ < Ce, from which 
we conclude that T vanishes whenever all Tj vanish. It follows from a standard argument 
(for example, Lemma 3.9 of |28 1) that T is a linear combination of the Tj. 

The argument described in the previous paragraph motivates us to find conditions on 
a test function <p that ensure we can cut if off outside a small neighborhood of a point q 
while keeping the norm \<p q \„, of the cutoff small. In order to proceed we will need some 
notation for a neighborhood base of q. If q is a non-junction point then it lies in a single 
copy of X in the cellular structure, and within this copy there is a unique word w such that 
F W (X) = q. The cells containing q are then of the form {/, = F [ W ],(X). For junction points 
the situation is different, as q can be the intersection point of several copies of X, or can be 
in a single copy but be given by F W .(X) = q for a finite number of words w\,...,Wj. We 
will not distinguish between these possibilities but will instead make the convention that 
the distinct words determining q may be used to distinguish copies of X if necessary. With 
this assumed, let Uij = F[ W ^{X), and U, = U 7 t/,j. 

Fix q and let Gy denote the Dirichlet Green's operator on the cell Uij, omitting the j 
index if q is a non-junction point. If <p e £>(Q) we can then decompose <p on Uij into 

,h- ( -//, ;C 4-r;,,V4 

where Hij<f> is the (unique) harmonic function on £/y whose values on dUjj coincide with 
those of . By induction we obtain 

m-l 
1=0 

and write h\ . = G\ fli iA'd. 

hj hj ' J 

Lemma 6.1. In the decomposition ( 16.11 ) we have at each x € C7y and p e <9t/,'j that 
A k h[.(x) = d„A k h[ J (p) = 0ifk> I, while for k< I, 

|A%|<c(*,/)rK ]( ^ ]( ||HyAV|| L . (l/y) 

\d n A%(p)\ < C (M)rf4X4IK^L Wj) 

Proof. The conclusions for the cases k > I are immediate from the fact that A l h\ . is har- 
monic, while the remaining estimates are derived from the fact that the Laplacian scales by 
r w fi w on a cell F H .(X) while the normal derivative scales by r w . □ 

Our purpose in making the above definitions is that estimates on the functions h l { . are 
precisely what is needed to ensure we can cut off a smooth function in the manner previ- 
ously described. 



DISTRIBUTION THEORY ON P.C.F. FRACTALS 



17 



/ \m-l 

Theorem 6.2. If ip is such that A' n cf>(q) - and WHijA'tpW^^^ = oUtw/l/V/l) f or 
< I < m - 1 as i — > oo, then for all e > there is <p q such that \(f> q \ m < e and (f> — <p q is 
supported away from q. 

Proof. We begin by constructing a neighborhood of U,j by adjoining cells at each of the 
points p e dU,j. At each p we require finitely many such cells, and we choose them so 
as to intersect Uu only at p. It will also be convenient to assume that these cells have 
comparable scale to the Uu, in the sense that they have the form F„(X) for some word 
with length \W\ < i + z'o for some constant i'o. Let K be one of the cells adjoined at a point 
p, and let n p be the number of cells adjoined at p. Using Theorem l2.3l we define a smooth 

function f K on K with jet A k f K (p) = A k <p(p) and d%A k f K (p) = -(l/n p )d" IJ A k (f>(p), and 
with vanishing jets at the other boundary points of K. Having done this for the set 7C of 
adjoined cells we see from the matching conditions for the Laplacian that 



?0) = 



(p(x) for x e Ui 

otherwise 



We must estimate |0J m . There is an easy estimate for A k (f> for k < m from Lemma loTTl 



defines a test function with the property that - <p q — on {/, 
must estimate \tp q \ m . There is an i 

m-l 

|Afy| < 2 \& k h l u \ + \G'lj k A m (j)\ 



m- 1 



<£c(fcj)r[;4^ 

l=k 



< 

l=k 

~(-jn—k ,jn—k\ re n\ 

where we used A"'(p(q) — to obtain that A m (f)(q) — o(l) when z — > oo. As a result we have 
good control of \<p q \ m on Ujj. 

A similar calculation allows us to estimate the size of the normal derivative \d„A k <f>(p)\ 
at any of the points p where pieces fx are attached. We compute 

m-l 

\d n A k </>(p)\ < Y) d " A %\ + \dnG™- k A m <f>\ 



m-l 



< J] c(k, /»iV/^ ll'/' ^IL „ , + c(k, m)^4-^4||A» 

l=k 
m 

<J]c(k,t)o(r?^p.?- k ) 

l=k 

Fix K e % and examine f^. By assumption = F%(X), so by ( 12.3b with the fixed 
number of jet terms m we know 

m m-l 

||A*/rlU < C(k)(Y 4- k 4- k \A k '<f>(p)\ + J] rl +1 - k 4- k \d^A k '<f>( P )\) + * (6-4) 

k'=0 k'=Q 
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provided < k < m. The terms involving |A <p(p)\ may be replaced by the estimate ( 16.21 ). 
For the terms involving normal derivatives we use that d„A k h k j(p) = (l/n p )d„A k <p and 
(16.3I>. The result is 



k'=o k>=a 

However, w and [vv,],- have comparable length and are adjacent, so they differ only in 
the final z'o letters and therefore the ratios r^r^^ and fiftf/.^^ are bounded by constants 
depending only on io and the harmonic structure and measure. It follows that 

lIAV/dL < C(m,k,r,p.)o(r [u!j]i _ l fi [H . j]i _ l )'"~ k 
and combining this estimate for each K e'K with ( 16.21 i proves that 

I|A*0 ? ||l« = o(r [H)j .] w jU [w . ]._,)"*" as i -> oo 
for < k < m. In particular we can make \<p Cj \ m < e by making i sufficiently large. □ 



Theorem l6.2l suggests that the natural candidates for the distributions supported at q are 
appropriately scaled limits of the maps <p m> Hi jA'<p as i — > oo. The question of how to take 
such limits has been considered by a number of authors 1 16, 38, 25, 26, 1 ], and is generally 
quite complicated. At the heart of this complexity is the fact that the local behavior of 
smooth functions in a neighborhood of a point q depends strongly (in fact almost entirely) 
on the point q rather than the function itself. This property - often called "geography is 
destiny" - contrasts sharply with the Euclidean situation where neighborhoods of points 
are analytically indistinguishable. Its immediate implication for the structure of distribu- 
tions with point support is that the nature of these distributions must depend on the point 
in question. In order of increasing complexity we consider three cases: junction points, 
periodic points and a class of measure-theoretically generic points. 

Junction Points. As before, the junction point q is q — F W .(X) for words w\ , . . . wj, each 
of which terminates with an infinite repetition of a single letter. The distributions corre- 
sponding to approaching q through the sequence [w,], may be understood by examining 
the eigenstructure of the harmonic extension matrices A, , the definition of which appeared 
in the Harmonic Functions part of Section|2] 

For notational convenience we temporarily fix one contraction F, let A be the corre- 
sponding harmonic extension matrix, and suppose q is C\F'(X). Let r and /u be the re- 
sistance and measure scalings of F, and ji, ■ ■ - ,y n be the eigenvalues of A, ordered by 
decreasing absolute value, with eigenspaces E\,...,E„. Of course j\ — 1 and E\ is the 
constant functions. Let H{u be the harmonic function on F'{X) that equals u on dF'(X), 
and P s be the projection onto E s . In what follows, G is the Dirichlet Green's operator on 
X and G, is the same on F'(X). 

Deflntion 6.3. Inductively define derivatives d s and differentials D k , k > 1 at the point q 
by setting D°u = u(q), and for each s such that (rfi) k < y s < (r/i) k ~ l 

d s u = lim y-JP s Hi(u - GD k l Au) (6.6) 



DISTRIBUTION THEORY ON P.C.F. FRACTALS 



19 



if these limits exist. Note that d s always exists for harmonic functions as the sequence is 
constant in this case. Provided the necessary d s u exist we then let 

D k u - h + GD kl Au (6.7) 

where h is the unique harmonic function on X with d s h = d s u for those s with y s > (rp) k 
and d s h = for all other s. We will also make use of D k u, where D° = m(0) and 

D k u = h + GD k ~ l Au 

where h is harmonic on X with d s h = d s u for those s with y s > (rfj) k and d s h = for all 
other s. 

Lemma 6.4. For u e dom(A 4 ) and each s with y s > (rfj) k the derivative d s u exists, and 

k 

Kw|<C(fe)^||A ; M |U. (6.8) 

;=o 

The differential satisfies 

\\u - D k u\\ mx)) < C(k)i k (r^) ki \\A k U \U (6.9) 

and if we further suppose that A k u 6 dom(£') then 

\\u - D k u\\ LX(F , (X)) < C{k)( m fS l2 E x l\A k u). (6.10) 

Proof. The proof is inductive. When k — there are no s with y s > 1 = {rp) Q , so the first 
statement is vacuous and J6.9I ) is immediate. Suppose both hold up to k - 1. 
Write u - GD k ~ x Au as H Q u + G(Au - D k ~ l Am), from which 

d s u = d s H u + limy^P s Hi(GAu - GD k ~ l Au) (6.1 1) 

i 

provided the latter limit exists. On the cell F'(X), 

G(Au - D k - l Au) = HiG(Au - D k ~ l Au) + G,(Am - D k ~ l Au) 

thus 

H M G{Au - D k ~ l Au) = AHiG(Au - D k ~ l Au) + H M Gi{Au - D k ~ l Au). 

In particular, if we project onto the eigenspace E s then the action of A is multiplication by 
y s . Scaling implies G;(Aw - D k ~ l Au) is bounded by 

\Gi(Au - D^'Am)|< C(rfxy\\Au - D^ 1 Au\\^ (F , (X)) < CC(k - ^'(r/^HA^IU (6.12) 
and the action of and P s can only improve this estimate, so 

y- {M) \p s H M G(Au - D k ~ l Au) - y s P s HiG(Au - D* -1 A«)| 
<y- (M) \Gi(Au-D k - x Au)\ 

^CCCt-l^f-^VllA^IU (6.13) 
v y s ' 

This shows {yJ'P s HjG(Au - D k ~ l Au)} is Cauchy when y s > (r/i) k , and that its limit is 
bounded by C(A;)||A' : m|| 00 . It follows from ( 16.111 ) that d s u exists for these values of s, and 
since \d s Hou\ < we also obtain ( 16.8b . 
Summing the tail of ( |6.13t establishes that 

k k 

\d s u - y'JPsHiiu - GD k - x Au)\ < C(k)i^—)\\A k u\\ 00 . 
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Now let h be the unique harmonic function with d s h = d s u for those s with y s > (rp) k and 
d s h = otherwise. Since yJ'P s Hjh = d s h is a constant sequence, we find 

\P s Hj(u -h- GD k - l Au)\ < C(k)(rp) ik \\A k u\U (6.14) 

for those s with y s > (rp) k . Recalling D k u — h + GD kl Au from (16.7b write 



(u - D k u)\ F . = Hi(u - D k u) + G,(A(m - D k u)) 

Ha.a , n { a .i r>*-l. 



= //,(« - /i - GD Aw) + G,(Am - D k ^Au). (6.15) 

We have estimated G,(Am - D k ~ l Au) in ( 16.121 i and the terms P s Hi(u - h - GD k l Au) for 
y s > (rju)* in ( 16.141 ). What remains are the terms P s Hi(u - h - GD kl Au) for y s < (rpj) k . 
Each of these is obtained as a sum, with 

i-i 

\P s Hi(u -h- GD k ~ l Au)\ = V y'^PsHjGjiu -h- GD k ~ l Au) 



./=» 

i-1 



./=<> 

i-1 

< CC(jfc- 1)(^)*||A*«|U V Mj^rf 

< C{k)i k (r^) ik \\A k u\U (6.16) 

because y s < (rfi) k . This proves J6.9t for and completes the induction. 

The proof of ( 16.10b uses essentially the same inductive argument with D replacing D 
and the estimate from ( 16.10b replacing that from d6.91 l throughout. Note that in ( 16.13b 
we can have y s > (r/u) k because there is an additional factor of r' /2 so the series still 
converges geometrically. Also, in ( 16.161 ) the working is simplified because for D we have 
these y s < (r/u) k and the r'^ 2 term is bounded, so the convergence is geometric here also. 
This allows us to remove the polynomial term in i. The base case k — is true because of 
the Holder estimate (12.21). □ 



The map d s takes a smooth function to the eigenspace E s , We now fix orthonormal 
bases for each of the E s , and refer to the co-ordinates of d s with respect to the basis for E s 
as the components of d s ; these components have values in C. 

Corollary 6.5. Each component d sv of a d s for which (rfj) k < y s < (rp) k ~ l is a distribution 
supported at q and of order at most k. Ify s < (rp) k ~ l then its order is equal to k, and it is 
otherwise of order either k — 1 or k. If d sv is a component that is a distribution of order k, 
then A d s>v defined by = d StV A l <p is also supported at q and has order k + I. 

Proof. It is apparent from the definition that d s is linear on D(Q.) and that d s (f) — if 
<p e £)(Q) is identically zero in a neighborhood of q, so it follows from ( 16.8b that the 
components of d s are distributions of order at most k and are supported at q. 

Suppose y s < (rp) k ~ l and let u s v denote the harmonic function determined by the eigen- 
vector corresponding to d s r . Then the values of u sy are 0{y s ) = o{rp) k - 1 and A l u = for 
/ < 1, so Theorem 16.21 implies that for any e > there is a function if/ equal to m SjV in a 
neighborhood of q but with \i[s\k-\ < £■ Since d SlV u SiV = 1 and d s , v u SiV = d s<v if/ because of the 
support condition, it cannot be that d s , v is order k - 1 or less, so it has order k. 

In the case y s = {rp) k ~ l < (rp) k ~ 2 the above argument says that d s , v has order at least 
k - 1. Both of the values k — 1 and k occur in examples. For instance, when k = 1, the 
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derivative d\u = u(q) corresponding to the constant harmonic functions has order = k— 1. 
A case where there is a d s of this type with order k occurs on the Sierpinski Gasket, see 
Example 16.71 below. This shows that scaling alone cannot identify the order of d s when 
y s = (rfif- 1 . 

The statement regarding A'd s is immediate. □ 

We now return to using the index j to distinguish the words Wj for which x = F Wj (X), 
and accordingly denote by d{ the derivative d s corresponding to the approach through cells 

Theorem 6.6. Let T be a distribution of order k supported at the junction point q, where 
q — F Wj (X), j — 1, . . . , n. The word Wj terminates with infinite repetition of a letter which, 
by a suitable relabeling we assume is j. Then T is a finite linear combination of the 
distributions A l d^ v , for which y, > (r^i,)* - '. The linear combination runs over all such s, 
all basis elements vfor E s , and all cells j — 1 , . . . , n that meet at q. 

Proof. Suppose that (p e D(Q) has the property that A l d J StV <t> = for all (rjp.f) k ~ l < y s . It 
follows from Definition ^. 3l that D k <p = and more generally that D k ~ l A l <p - for all I < k. 

However, the harmonic part of HjjA'cp = H,jA l ((p - D k <p) on the cell C/y of scale i 
corresponding to the word w j is bounded by the maximum over the boundary vertices of 
this cell, so from (16.10b : 

IK A VL (£/y) = oirjVjf-* = o{r [wmWjl ) k -' 

for < I < k - 1. We also have that A k <p(q) = because A'^(^) = A l d[ v <p = 0, so 
Theorem 16.21 shows that for any e > there is iff E £)(Q) that is equal to <p - D k <p in a 
neighborhood of q and with 1^% < e. 

Using the support condition and the fact that T has order k yields 

T(p = T\p < M\i//\ k < Me 

for some fixed M depending only on T, and all e > 0. Thus Tip = 0, and we have 
shown that the kernel of T contains the intersection of the kernels of the distributions 
described. By a standard result (e.g. Lemma 3.9 of [28 1), T is a linear combination of 
these distributions. □ 

Remark 6.7. Since d\ corresponds to the eigenspace of constants, the distributions d[ A 1 
are independent of j and are simply powers of the Laplacian applied to the Dirac mass 
at x. It should also be noted that for each j the distribution dl corresponds to the largest 
eigenvalue less than 1, so gives the normal derivative at x when approaching through the 
cells F[ W .]i, i — > oo. As a result 2 j d J 2 u = 0, and not all of these distributions need appear 
in 7". 

It should also be noted that the linear combination in Theorem l6 . 6l mav include distribu- 
tions of the form A l d J StV having y s = (rjfij) k ~ l , and that it is possible for these to be of order 
k + 1 . If this were to occur then we would have a non-trivial linear combination of these 
(k + l)-order distributions such that the linear combination is of order only k. We do not 
know of an example in which this occurs, but cannot eliminate it as a possibility because 
our arguments rely on scaling information. 

Example 6.8. The canonical example of a p.c.f. self-similar fractal of the type we are 
describing is the Sierpinski Gasket SG with its usual symmetric harmonic structure (see 
ll36l for details of all results described below). In this case r = 3/5 and yU = 1/3, so the 
Laplacian scales by rp = 1/5. Each of the harmonic extension matrices A, has eigenvalues 
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1, 3/5 and 1/5 with one-dimensional eigenspaces. The corresponding derivatives at q are 
ofjM = u(q) which is point evaluation, d J 2 u(q) = d J N u(q) which is the normal derivative at q 
from the cell corresponding to j, and d J 3 u(q) = d J T u(q) which is the tangential derivative of 
u at q from this cell. 

There are two cells meeting at the junction point q. Without loss of generality we 
suppose they are indexed by j = 0, 1 . The two corresponding normal derivatives d J N u(q), 
j — 0, 1 satisfy the single linear relation that they sum to zero, and by Corollarv l6.5l thev are 
of order 1. The two tangential derivatives d J T u(q), j — 0, 1 are independent, and it is known 
that they cannot be controlled by ||zt||oo + ||Au||oa (see ll36l . page 60). They are therefore of 
order 2. It is also possible to see in this example that any non-trivial linear combination 
of the d J T u(q) has order less than 2. Writing 5 q for the Dirac mass at q, we conclude from 
Theorem l6.6l that any distribution T of order k at a junction point of S G can be written as 
a linear combination of the form 

k k-\ k-2 

T = ^ ai A% + £ b,tid a N 6 q + c O' A '4V (6.17) 

/=0 1=0 1=0 j=0,\ 

This example also illustrates the issue described in the proof of Corollary 16.51 namely 
that there can be a d s with y s = (r(X) and yet d s is order k. In this case we have dj = dj, 
with 73 = 1/5 = rfi, so k — 2, and di is of order 2. 

Periodic and Eventually Periodic Points. Periodic points are those x = F W (X) for which 
w is a periodic word, meaning that w is composed of an infinite repetition of a fixed finite 
word v. Eventually periodic points are those for which the word w is periodic after some 
finite number of letters. For these points there is a theory similar to that used for junction 
points; we do not have to consider derivatives corresponding to multiple cells, but instead 
of looking at the eigenstructure of a matrix A, we must examine that of A v , which is a finite 
composition of the A, matrices. If y s is an eigenvalue of A,, with eigenspace E s , then we 
can define the derivative d s as we did for junction points. It is easy to see that the analogues 
of Lemma l6.4l Corollarv l6.5l and Theorem l6.6l all hold, simply by changing the notation to 
refer to the infinitely repeated matrix being A v , the eigenvalues y s being those of A v , and 
the Laplacian scaling factor to be r v fi v instead of rjfij. 

Generic Points. We now consider a non-junction point x = F W (X), where w = w\W2 ■ ■ ■ 
is an infinite word. The behavior of harmonic functions on the cell F[ W ] n (X) can be under- 
stood by considering the product A[ H ,j n = n?=i A Wj . We need to understand their scaling 
properties, for which we use the following approach from 11381 . Define for each unit vector 
a the corresponding Lyupunov exponent 

logy(a) = lim j log ||A Wj ar|| (6.18) 

if the limit exists. In this definition we may take || • || to be any norm on the #Vo-dimensional 
space containing a; all such norms are equivalent, so y is unaffected by this choice. 

Let us suppose that these limits exist at x. It is readily seen that y(a) + y{a') implies 
a and a' are orthogonal, so there are at most #Vb distinct values y\ > yi - ■ ■ that occur. 
Corresponding to these is a direct sum decomposition E\ © E2 ffi • • ■ with the property that 

writing a = a\ + aq. H we have y(a) = y s if and only if a\ — ■ ■ ■ — a s -\ - and a s + 0. 

Since the constant functions are harmonic we actually know that y\ = 1 and E\ is spanned 
by (1,1,..., 1). We let P s be the orthogonal projection onto E s . 

The subspaces E s provide the natural decomposition of harmonic functions into their 
scaling components at x. However we cannot expect to directly mimic Definition 16.31 
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because the estimate ( 16.181 ) does not imply the existence of a renormalized limit of the 
form 

lim y^P s Hi(u - GD kl Au) (6.19) 

i— >oo 

Indeed it is easy to see that ( 16.181 ) does not even imply that A^a is 0(y(a)) 1 . 

A natural way to proceed was introduced in [16, 38| and further treated in [26]. Let 
u be the function we are considering, and H/u be the harmonic function on F[ W ] t (X) with 
boundary values equal to u on dF[ W ] t (X) as usual. If we assume that the harmonic scaling 
matrices Aj are all invertible we can unravel the scaling structure for harmonic functions 
at x by applying the inverse of A[ W ], to H,u. For later use we record an elementary result 
about the scaling of the adjoint of Aj~Jp 

Lemma 6.9. If a e E s then lim i log||(A^)*a;|| = - logy s . 
Proof. Writing (•, •) for the usual inner product, 

» (A ^ g|1 = s 7 ikii 

\(A lw] ,a',(A-iya)\ 
= sup ■ 



\\A [w] ,a'\ 
\(a',a)\ 



= sup ■ 

a' \\A [w]j a'\\ 

Since the logarithm is monotone, this implies 

\(a',a)\ I \\A [w]i a'\ 



1 I, , . n (I m ,a)\ 1 

7 iog||(A^lh-p(7 log W-7 log - 



Of' 



however we know that the second term inside the supremum converges to -y(a'), whereas 
the first converges to zero provided (a', a) + 0. The latter condition and a e E s requires 
that a' have a non-zero component in E s , from which we deduce y(a') > y s , with equality 
provided P,a' = for each t < s. Combining these observations it is easy to see that 
for each i the supremum is between -y s - S and -y s for a constant c independent of i. It 
follows that the limit in the statement of the lemma exists and has the asserted value. □ 

In order to account for the scaling behavior of the Laplacian, we set 

logy6 lv . = lim — log r Mi yu Mi (6.20) 

provided the limit exists. 

Deflntion 6.10. Assume that x = F W (X) is a point at which the limits in d6.18t and ( 16.201 ) 
exist, and that all Aj are invertible. Inductively define derivatives d s and differentials D k , 
k> 1, at the point x by setting D°u = u(q), and for each s such that fi\ v <y s < /3„7 l 

d s u = lim P s A7*Hi(u - GD k l Au) (6.21) 

if these limits exist. Note that d s always exists for harmonic functions because the sequence 
is constant in this case. Provided the necessary d s u exist we then let 

D k u - h + GD k l Au (6.22) 

where h is the unique harmonic function on X with d s h = d s u for those s with y s > fi k w and 
d s h = for all other s. We will also make use of D k u, where D° = m(0) and 

D k u = h + GD k - l Au 
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with h harmonic on X with d s h = d s u for those s with y s > j3 k w and d s h = for all other s. 

Observe that this generalizes Definition 16.31 because if x = F W (X) is a junction point 
then w ends with infinite repetition of a single letter j, the Lyapunov exponents are the 
eigenvalues of Aj, and the action of Aj"^ on the eigenspace E s is just multiplication by y~'. 

The following result may be seen as a generalization of Theorem 1 of [38], see also 
Theorems 5 and 6 of [26 1 . It is proved by essentially the same method as Lemma l6~4l At 
several points in the proof we use the observation that for a positive sequence a, satisfying 
limi" 1 logo, = logo and a value e > there is a constant C(e) so C(e) _1 e _e! a ! < aj < 
C{e)e ei a l . 

Lemma 6.11. Assume that all Aj are invertible, and that x — F W (X) is a point at which 
the limits in ( 16.181 1 and d6.20| l exist. For u e dom(A*) and each s such that y s > /?*., the 
derivative d s exists, and 

k 

\d s u\<C(k)Yj\\^u\U (6.23) 

1=0 

For all sufficiently small e > 0, the differential satisfies 

\\u - £>*m|Uf [u , ( x)) ^ C ^ IIA^IU (6.24) 

If in addition we assume that A k u e dom(£') then 

\\u - D k u\\ L ^ Fm{x)) < C(k, e)r\ll0yE X! \A k u). (6.25) 

Proof. The proof is inductive. When k — there are no s with y s > 1 = so the first 
statement is vacuous and ( 16.241 i is immediate. Suppose both hold up to k - 1. 
Write u - GD k l Au as H Q u + G(Au - D* -1 Ak), so 

d s u = d s H u + lim P s A[^Hi(GAu - GD k l Au) (6.26) 

provided the latter limit exists. Writing G, for the Dirichlet Green's operator on the cell 
F[ W ]i(X), we have on that cell 

G(Au - D k ~ l Au) = HiG(Au - D k ~ x Au) + G,(Am - D k ~ x Au) 

from which 

H M G{Au - D k ~ x Au) = A WM HiG(Au - D* _1 Ak) + H M Gi(Au - D fc_1 A«), 
therefore 

A^ ]M H i+ i(GAu - GD k l Au) - A^H^GAu - GD k ~ l Au) 

= A-^H^GiiAu - D k - l Au), (6.27) 

and by substitution into ( I6.261 l. 

oo 

d s u = d s H Q u + ^ PsA^H^GiiAu - D* -1 Aw) (6.28) 
o 

provided that the series converges. 
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Since G,- inverts the Laplacian with Dirichlet boundary conditions on F[ W ] ( (X), we have 
for any sufficiently small e > the bound 

|G;(Aw - D* -1 Ak)| < Cr [H . mw]i \\Au - D^Aw^^ ) 

<C(*-l,e)$*< e ^||A^|| M (6.29) 

because of the inductive hypothesis ( 16.241 ) and the Laplacian scaling estimate ( 16.201 ). This 
is also applicable to //,+iG /(Am - D k ~ l Aw) by the maximum principle. Using Lemma |6T9l to 
estimate the size of \\(A^ )*P s a\\, it follows that for any sufficiently small e > 0, and any 
vector a, 

^P^HMGiiAu - D^Au), a)\ = \{H i+l GiAu - D k ~ l Au), (A^ ]( )*J>,a>| 

<C(fc-l,6)^ r r e (3e/4, 1|A«|L, 

This and the assumption y s > f3 k v imply that if e > was chosen small enough then 
the series in (16.281 ) converges, and is bounded by CIIA^mIIoo. The estimate (16.231 ) follows 
because d s Hoii is bounded by C||m||oo- 

Now u-D k u — u — h- GD Au = Hqu - h + G(Au - D 1 '' 1 Au), where h is the harmonic 
function with d s u = PJi for all s satisfying y s > /3 k K and P s h = otherwise. An expression 
for h can be obtained by summing ( 16.281 ) over these values of s. Comparing it to the 
expression 

A^Hiiu - D k u) ^Hou-h + Yj A^HmGAAu - D k ~ x Au) 

o 

from ( 16.27b . it is apparent that for those s with y s > fi k w we have 

oo 

PsA^Hiiu - D k u) = - J] P s A^H M Gi(Au - D^Au) 

i 

which we note satisfies for all \\a\\ < 1 and sufficiently small 6 > 

CO 

\{P s A^.Hi(u - D k u), a)\ < Y){H l+1 G,(Au - D^Au), (A^ h yP s a)\ 

i 

oo 

<Y,C(k-l,e)/3 k J,y;'e^%A k u\l 

i 

<C(jt-l,e)^y; i <! <3e/4) '||A*8|| oo . (6.30) 
For those s satisfying y s < p k w we have instead 

i 

PsA^Mu - D k u) = J] P s A^H M Gi(Au - D k ~ l Au). 
o 

and for all vectors a with ||ai| < 1, 

i 

{(PsA^Hiiu - D k u),a)\ < J] C(k - 1, e)/£V;'e (3e/4) i| A *"L 

o 

<C(^-l,6)^ r ;' e (3e/4)i ||A*«|| M . (6.31) 
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Equations ( 16.301 ) and ( 16. 31b give the same estimate for each P s A^Hj{u - D k u). Mapping 
forward again by A[ w ] t increases each term by a factor at most C(e)yJ s e <£ ^ 4) ', so summing 
over all s we finally have 

\H t (u - D k u)\ < C^"e ei ||A*M|| oo 
for some constant C - C(k, e). Now the restriction of (u - D k u) to F^ v y(X) is 
(u - D*w)| = Hi(u - D k u) + G,(A(m - D k u)) = H^u - D k u) + G t (Au - D k ~ x Au)) 

the second term of which is bounded by /j{jfe s ||A*H|| from ( 16.291 1, and the first term of 
which we have just estimated in the same way. This establishes ( 16.24b and completes the 
induction. 

The proof of (16.25b is the same, except that ( 16.251 is used in place of ( 16.241 > throughout. 
The validity of the estimate for k — is a consequence of the Holder estimate (12.2) . □ 

As previously, we fix orthonormal bases for the spaces E s and see that the components 
of d s are distributions. 

Corollary 6.12. Suppose that x satisfies the assumptions of Lemma \6. 1 1 \ and /3 k . < y s . Any 
component d s , v of the derivative d s is a distribution of order at most k supported at x. If 
also y s < then d StV has order equal to k. If d SyV has order k then defining A l d SyV by 
Ad s>v (p — d s>v A l <j> yields a distribution supported at x and of order k + I. 

Proof. Linearity of d s , v is immediate from Definition 16. 101 so it is a distribution of order 
at most k by (16.23b . Again using Definition ^. lOl it is apparent that d StV (p — if <p £ !D(£l) 
vanishes in a neighborhood of x, so d SiV is supported at x. 

To see that d s , v has order at least k, consider the harmonic function h with boundary 
values equal to the unit vector in the v direction in E s . Then HJi = A^Hgh, so the 
sequence in (16. 2U is constant equal to Hoh, and d s v h = 1 . Now for e > so small that 
y s e 3e < p k y X we have 

H^IL < C(e) r ^ ri < Cie^-^e- 2 " < C{e){r Mi p [w]i ) k - 1 = o{r Mi p [w] t l 

and of course A l h = for all I > 0, so Theorem l6. 21 applies with m — k - 1, and there is a 
test function <p such that (f> = h in a neighborhood of x and \<p\k-\ is as small as we desire. 
Since d s , v h = d StV <p by the support condition, d s v cannot be of order k - 1 or less. The final 
statement of the lemma is obvious. □ 

Theorem 6.13. Suppose that all of the matrices Aj are invertible, and that x — F W (X) is a 
point at which the limits in ( 16.18b and ( 16.201 ) exist. Then all distributions of order at most 
k at x are linear combinations of the distributions A d sy , with y s > p£ . 

Proof. As in the proof of Theorem 16.61 it suffices to show that T vanishes whenever the 
distributions A l d s , v , with y s > fi k ~ l vanish. 

Suppose <p e D{Q) satisfies A l d StV <f> = for those y s > j3 k ~ l . Then the differential 
D k (p (which exists by Lemma 16.1 lb must be zero, as must D k ^A l (f) for each < I < k. 
From ( 16.25b we then see that for all sufficiently small e > 0, 

II^U Mi (*)) £ ll A ^IU M ,<*)) * C(k,e)rlilpf-W'\A k u) 

< C{k, e)r^(r Mj /i Mj )*-'e M £ 1/2 (A i «) 
= o(r Mi fx [w]i f~ l . 
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Applying Theorem l6.2l we find that for any 6 > there is i// equal to in a neighborhood 
of x and such that \\fz\k < 5. In particular, since T is order k and supported at x, there is M 
independent of (f> such that 



In concluding this section it seems appropriate to say something about the set of points 
x satisfying the conditions in Definition ^. 101 The set at which the limit f} w exists has full 
//-measure by the law of large numbers, and in fact 



at /i-a.e. point. The set on which the Lyupunov exponents exist may be treated by the 
theory of random matrices introduced by Furstenberg and Kesten 0. In particular, it is 
possible to make certain assumptions on the matrices A, that guarantee that this set is also 
of full //-measure. This topic is discussed quite thoroughly in the paper [26] of Pelander 
and Teplyaev, so we will not cover it here. One consequence of their work, however, is 
that there are conditions that imply the spaces E s are independent of the choice of point x. 
For example, if the semigroup generated by the A, is strongly irreducible and contracting 
then there is a single vector a\ such that at //-almost every x, the space E\ is spanned by 
a\ and has scaling j\ . If the same strong irreducibility and contraction holds after taking 
the quotient to remove E\ then E% is also one-dimensional and independent of x on a full 
measure set. For a fractal where the irreducibility and contraction properties are true for the 
semigroup generated by the A, on each of the subspaces found by removing E\,E2,.., E s -i 
in turn, we could conclude that all of the distributions of the form d s are independent of x 
on a set of full //-measure. Hence in this situation any distribution of order m with point 
support in a fixed set of full //-measure would be a finite linear combination of distributions 
A'd s for suitable values of /, where the d s are independent of x. This generic behavior is 
very different from that seen at junction points and eventually periodic points, where the 
structure of point-supported distributions can vary substantially from point to point. 



In this section we give a theory of distributions on finite products of post-critically finite 
self-similar fractals, using the analytic theory for such products developed in [34|. This 
gives genuinely new examples, because products of p.c.f. self-similar sets are not usually 
themselves p.c.f. Since there is no essential difference between a products = X' XX" with 
two factors and a general finite product, we state our results only for the two factor case. 

Following the notational conventions of l34l . points are x — (x', x"), functions on X are 
called u or /, on X' they are u' or /' , while on X" they are u" or /" . The energies on X' and 
X" are E' and E" and the Laplacians are A' and A". They come from a regular harmonic 
structure as in Section [2] and have energy and measure scaling factors r', //', r" and //". 
The corresponding Laplacians A' and A" are defined componentwise, so u e dom(A') with 
A'u — f if u and / are continuous on X and have the property that for each fixed x" e X" 
we have A'u(-, x") = /(•, x"). A similar definition is used for A". By Lemma 1 1.2 of 041 . 
A' and A" commute on dom(A') n dom(A"). 

Deflntion 7.1. A function u on X is smooth if for all j, k eN, (A') J (A") k u is a continuous 
function on X. The definition extends to a finite union of cells in the obvious manner, and 
u is smooth on a domain in X if it is smooth on every finite union of cells in the domain. 



\T<p\ = \Tif,\ < M\\fi\ < 6 



and thus Tcp = 0. 



□ 




j 



7. Distributions on products 
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We define the test functions on a domain D. to be the smooth functions of compact 
support with the usual topology and the corresponding seminorms 

|0U = sup{|(A'y'(A")VW| : x e Q, j + k < m"). 

The distributions form the dual space with weak-star topology. A distribution T has order 
m if on any compact K there is M — M(K) so that \T<p\ < M\<f>\ m for all test functions 
supported on K 

The goal of this section is to provide conditions under which analogues of our main 
results for distributions on p.c.f. fractals are also valid on products of these fractals. In 
order to avoid duplicating a great deal of work, we only give details of the proofs where 
they differ significantly from those for the case of a single p.c.f. fractal. In particular it is 
fairly easy to verify that all of the results of Section[3](except Corollarv l3.61 l. SectionH] and 
Section [5] prior to Theorem 15. 61 depend only on the partitioning property of Theorem l2.7l 
and the estimate (12.5b (either directly or through Lemma [3~3l as well as the fact that for 
compact Q there is an orthonormal basis of l? consisting of eigenfunctions. The latter is 
obviously true for the product X' x X" because it is true for the factors, so the original 
proofs transfer to the product setting once we know the partitioning property and the cor- 
responding estimate for product spaces. These are proved in Theorem l7.7l and (17.6b below, 
so all of the aforementioned results are also true for products of p.c.f. fractals with regular 
harmonic structure, and connected fractafolds with restricted cellular structure based on 
such products. 

Only small changes are needed to obtain analogues of the remaining results from Sec- 
tion [5] The proof of Theorem 15.61 required that on any cell there was a positive smooth 
function equal to 1 on the cell and vanishing outside a specified neighborhood: such a 
function may be obtained on the product space as a product of functions of this type on the 
factors, so the theorem is true for products in which each factor has the estimate ( 12.4b for 
the heat kernel. The other results are used to prove the structure theorem (Theorem |5.9l ). 
Of these, Lemma l5?7l remains true with the same proof if it is modified to say that A' maps 
D{K) to itself with image orthogonal to those (f> having A'(f> = and there is G' K such that 
-A'G' K (A'(p) = A'0; there is a corresponding result for A". A version of Theorem l5.8l is 
then true with G'^G' K replacing Gk throughout. The original proof shows that for m > 1, 
G'^G' K takes a distribution of order m to one of order at most m — 1 . To show that G'^G' K 
takes a distribution of order zero to a continuous function it suffices to approximate the 
corresponding measure v by a sequence of linear combinations of product measures. Ap- 
plying G'^G' K to a product measure gives a continuous function by the original proof of 
Theorem 15.81 so applying it to the sequence gives a uniformly convergent sequence of 
continuous functions whose limit represents the distribution G'^G' K v. The proof of Theo- 
rem l5.9l needs no further changes. 

At this point we have essentially all of the results of Sections [3] 111 and|5]in the product 
setting (the only exception is Corollarv l3.6b . In addition there are some things that can be 
said about distributions with point support that generalize the results of Section[6] We will 
return to these after giving the details of the partitioning argument, because some aspects 
of the procedure for cutting off a smooth function will be important for the proofs. 

Partitioning on products. We prove analogues of the partitioning property in Theorem l2.7l 
and the estimate ( 12.5b in the product setting. As in the single variable case, the proof re- 
lies on a cell-by-cell construction of a smooth function, for which the following matching 
condition is essential. Note that a cell in X is a product of cells from X' and X", so has the 
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form K = F' W ,(X') x F'^,,(X"), where w' and w" are finite words. Its boundary consists of 
faces {<?;} x F»,(X") and F'JX') x [q'j] for q'. e V' and q>! e V^'. 

Lemma 7.2. Suppose the cells K\,...,Ki all contain the face L — {q'\ x F",„(X"), and 
that the union Uj'A'/ contains a neighborhood of every point in L except those of the form 
(q', F'^„q") with q" e Vo- If u j is smooth on Kj for each j, then the piecewise defined 
function u — Uj on Kj is smooth on U^K/ if and only if for each x" e X", both 

(a) The functions {A')\A") m Uj(q', x") are independent of j for each I and m, and 

(b) For each x", Y J j(d'„)j(A') l (A") m Uj(q' , x") — 0, where (d' n )j indicates the normal deriv- 
ative in the x' variable from within Kj. 

Proof. For fixed x", (b) is the necessary and sufficient matching condition in the first vari- 
able for (A')'(A")" ! m(-, x") to exist (as a function rather than a measure with atom at q'). 
Condition (a) is then equivalent to continuity of (A')'(A") m u. a 

Our construction uses an analogue of the Borel theorem from [27]. That result yields 
the existence of a smooth function with a prescribed jet at a junction point of a pcf fractal, 
whereas we need existence of a smooth function with prescribed smooth jet on the face of 
a cell in the product X. 

Theorem 7.3. Fix a face {q'}xX" and a neighborhood U c X' of q' . Given two sequences 
jp;(x")}" () and {o";(x")}" of functions that are smooth in x" , there is a smooth function u 
with support in U x X" such that for each x" e X", (A') k (A")"'u(q',x") = (A") m p k (x") 
andd' n (A'f(A") m u{q',x") = (A") m o- k (x"). 

Proof. The proof is almost the same as that for Theorem 4.3 of 1271 . Specifically we form 
the series 

u(x',x") = ^pi(x")g hmi (x') + o-i(x")fi Ml (x') (7.1) 
I 

where the functions and // „, are as defined in that proof, so they satisfy 
(A')W<?) = hi KiA'fg^q) = 

(A') k f L „,(q) = d' n {A') k f Uni {q) = Su 

and have supports in cells of scale m/ and «/ respectively. Convergence of the series d7.U 
is achieved by making an appropriate choice of m; and «/. In particular, it follows from 
the cited proof that if |p/(x")| < Ri and |cr/(x")| < Si for all x" , then one can choose 
mi and n; depending only on R/ and 5/ such that for each x" the series converges to a 
function that is smooth in x', supported in U x X", and has (A') k u(q, x") = pk(x") and 

d' n (A') k u(q,x") = cr k (x"). 

Now we require convergence not only of the series for u(x',x"), but also that for 
(A") m u(x' , x") for each m, so we must diagonalize. Set 

Ri = max max|(A")' Pi(x")\ 

0<l"<lx"eX"' 1 

Si = max max|(A")' o~i(x")\ 

0<l"<l.x"eX" 1 1 

which are finite by the assumed smoothness and the compactness of X" , and let m; and «/ 
be chosen as described above. For fixed x", all terms after the m-th in the partial sum 

L 

(A")'" J]pi(x")gi, mi (x') + o-i(x")f Uni (x') 
i 
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have coefficients bounded by Ri and S so the above reasoning implies that the partial sums 
converge to a function that is smooth in x', and has (A') k (A")'"u{q', x") = (A")'"p*( x ") an d 
d' n (A') k (A") m u(q\ x") = (A") m o- k (x") for all m. 

Finally, it will be useful later to have estimated the contribution of each term to the L°° 
norm of (A') k (A") m u. It is convenient to write w'(;«/) and w'(ni) for the words such that 
F w >( m ,)(X') is the support of gi^ and F' w , (n ) (X) is the support of /;_„,. Note that scaling then 
implies (see equations 4.4 and 4.5 of J27|) that 

\\{A')\A'T P i{x'')gi, mi ^)\\ < eft /)(C (m ,yu;, (m() )^||(A"rpKx")L 

||(A') fc (A"rcrK/')/ ; ,„ ; (^)|| < c(k, o(c ( „ ;) //;, ( „, ) ) / -v 1 ;„ ( „ )) ||(A") m C r ; (x")|L 

and in the construction in (271 it is noted that the contributions of terms with I > k may 
be made smaller than any prescribed e > 0, so taking e to be a small multiple of ||polU we 
obtain 

||(A')*(A' T «L < 2 eft o(C ( „,X'w)"ll( A ") m p'^")L 

1=0 

+ 2 c(k, 0(C ( „X. J|(ATVK*'')|L ( 7 -2) 

□ 

Remark 7.4. This result may be localized to any cell in X simply by rescaling the desired 
jet for the cell to obtain a corresponding jet on X, applying the theorem, and then compos- 
ing the resulting function with the inverse of the map to the cell. It may also be applied to 
a face in a finite union of cells, so that the face is of the form [q'} x (U . =V K") with each K" 
a cell in X", because \J J . =l K'! is compact. 

In order to make use of the preceding result we require a small lemma. 

Lemma 7.5. If u is smooth on X and q' e Vq then d' n u(q', x") is smooth with respect to x" 
and (A") l d' n u{q' , x") — d' n {A") l u(q' , x"). There is a bound 

\\d' n (A")' U (q', x")IU < C(||(A")'«|| M + ||A'(A")'«||J (7.3) 

Proof. For each x" and each scale m, let h,„(x', x") be the function that is piecewise har- 
monic at scale m in the x' variable and coincides with u on V' m x \x"}. Then h m (x', x") is 
smooth in x", because its values are obtained as uniform limits of linear combinations of 
the values from V' m x [x"). Moreover, the normal derivative d' n h m (q' , x") is a linear com- 
bination (with coefficients depending on m) of the differences {hmip'^x") - h m (p' 2 ,x")), 
where p' x and p' 2 are neighbors of q' at scale m. Thus d' n h,„(q' , x") is smooth in x" and 
(A")'d' n h m (q', x") = d'„(A")'h m (q', x"). 

For each fixed x" , we may express (A")'u(x', x") on a cell K' m of scale m containing q' 
as the sum of (A") l h m and an integral involving the Dirichlet Green kernel G' m for A' on 
K' m . Taking the normal derivative we obtain 

d'„(A") l u(q', x") = d'„(A")'h m (q', x") + j (A'(A")'«(/, x"))d' n G' m (q',y) dp'(y'). (7.4) 

However an easy scaling argument shows that d' n G' m (q' ,y') is bounded independent of m 
and y' , so the integral term is bounded by a constant multiple of HA' (A") z m|I ju'C-ST^), inde- 
pendent of m and x" . Since p'(K' m ) — > as w — > oo we conclude that (A") d' n h m {q', x") 
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converges to d' n (JS.")'u(q' , x") uniformly in x" for each I. Then ( 17.31 ) is obtained by us- 
ing (17.4b with m = 0. a 

We may use the preceding results to smoothly cut off a smooth function on a neighbor- 
hood of a cell. 

Theorem 7.6. Let u be smooth on a cell K = F' ,{X') x F",,(X"), and U D K be open. 
There is a function v such that v — u on K, v — onX\U and v is smooth on X. Moreover 
for each k, 

k m 

HtA'AA'TvL < C(k, U) £ £||(A'/(A'')«"|U/Q' < 7 ' 5 ) 

1=0 «=0 

Proof Let K' = F' w ,(X') and K" = F'^,(X"). Fix a face of K having the form {<?'} x K" 
and let pk(x") = (A') k u(q', x") and o~k(x") = d'„(A') k u(q', x"). The functions pi< are smooth 
in x" by the definition of smoothness of u, and the functions cr^. are smooth in x" by 
Lemma lTBl Now take a finite number of small cells K'j in X' with the following properties: 
the intersection K' n K'. = \q'} for all j, the intersection K'. C\ K'~. = {q'} for all j + j, the 
union K' U (UjK'j) contains a neighborhood of q' in X', and (K' U (U ; ^)) x K" c f/. Let 
the number of K'. be J, and apply Theorem 17. 3 1 to each K'. to obtain a smooth function uj 
that has jets pk(x") and (-l/7)cr^(x") at g ,/ and is supported in a neighborhood of q' that 
is strictly contained in K'j. By construction, the matching conditions of Lemma l7. 21 apply 
to the functions u onK and uj on K'j x /T", so the piecewise defined function is smooth on 
the union of these cells. 

Repeat the previous construction for each of the finite number of faces having the form 
\q'.} x K" . As these faces are disjoint we may choose the small cells in the construction so 
that those used for q' t do not intersect those for q'. for j + i. The result is a finite collection 
of cells K'j x K" c U and functions uj such that the piecewise function u on K and uj on 
K'j x K" is smooth on the union of the cells, and vanishes identically in a neighborhood of 
any boundary face of (K' U (UK'j)) x K" that has the form {p'} x K". We call this function 
v'. 

Having treated the vertical faces {q^} x K", we then treat the horizontal faces (K' U 
(UK'j)) x \q"} of the new function v' in the same manner. All of the results we needed were 
valid on faces of finite unions of cells, so the same proof allows us to piecewise extend to a 
smooth function v on a larger finite union of cells, which we call L, but with the additional 
condition that v vanishes identically in a neighborhood of each horizontal face of L. Then 
Let/ and v vanishes in a neighborhood of all faces of the boundary of L, so Lemma 17721 
ensures that extending v to be identically zero outside L gives a smooth function on X. By 
construction, v = u on K. 

For the estimate ( 17.5b we note that 

l|(A")W)L =s ||(A')*(A'THU*) 

by definition, while rescaling d7.31 > to the cell K implies that 

||(A")'W)L =S c((r;0i|(A')*(A''r«|U*) + ^||(A') t+1 (A"r«IL„ m ). 
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Substituting into (17.2b and using r' w ,^ n , < r' w , and fi' w , < 1 we have 

I|(A')*(A")1L * Z o(C (mi )^( m ,))"ll( A o'(A"r«IL 



(K) 



<C(^f/)2l|(A')'(A"r«|L w - 

/=0 

This type of estimate deals with all of the vertical faces, and an analogous argument is 
valid for the horizontal faces, so ( 17.5b holds. □ 

Theorem 7.7. If it is smooth on X and L)Q,j is an open cover ofX then there are constants 
C(k, m) and smooth functions uj such that uj is supported on Q.j, Yij 11 j — u > an d 

k m 

||(A')*(A") m ";IU < C(k, m) ^ J] ||(A')'(A")""IU (7.6) 

;=o n=0 

Proof. The open cover is finite, say because X is compact. Moreover we may par- 
tition X into a finite number of cells K/ such that each K[ is contained in some Slj. We 
proceed by induction on /, with the base case being that we apply Theorem l7.6l to u on K\ 
to obtain a smooth function vi with support in the open Q.j that contains K\, At the /-th 
step we apply Theorem l7.6l to u - £'f 1 v m on Ki to obtain a smooth function V; with support 
in the open Qj that contains Ki. Note that u - Yj\ v m vanishes on u'_ ^K m so once we have 
exhausted the cells we have Yui v i - u - By construction, each of the v; is smooth, supported 
on some Clj and satisfies ( I7.51 l. Setting uj to be the sum of those v/ that are supported on 
Q.j completes the proof. □ 

Distributions with point support on products. It is useful to begin with the observation 
that if T € D'(X') and T" £ 2D'(X") are distributions on the components of a product 
space X' x X" then there is a tensor distribution T x T" which is a distribution on the 
product. This is not entirely immediate, but follows readily from the structure theorem for 
the component spaces. Specifically, the fact that T is locally (-A') k f for a continuous / 
implies that for a <p e T>(X' x X") there are k and / such that 

A"T'(p(x',x") = A" \ f(x')(-A'f<fi(x',x")dfi'(x') 

JX' 

f(x')(-A') k A"<f>(x', x") dn'(x') 



L 



IX' 

= T'A"<p(x',x") 

where we used that A' and A" commute. In particular T'4> is smooth in the second variable, 
so T' x T"(f> = T"{T'<p) is well defined. Repeating the calculation with T" in place of A" 
ensures that T"(T'<f>) = T'(T"<f>), so the order in which the distributions are applied is not 
important. Linearity of T'xT" is immediate and it is easy to check the continuity condition 
that ensures it is a distribution on X' x X" . 

In the special case where T is supported at x' and T" is supported at x" it is apparent 
that T x T" is supported at (x', x"), so this construction and the results of Section|6]supply 
a large number of distributions with point support. In fact we can show that if x' and x" are 
either junction points or satisfy the conditions of Theorem l6.13l then the distributions with 
support at (x', x") are of this type. As in Section|6l the key is to show that if e D(X'xX") 
is annihilated by sufficiently large collection of tensor distributions at (x', x") and if e > 
is given, then it is possible to cut off on a small neighborhood of (x', x") such that the the 
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resulting function has (A'y(A")*0 < e for all j and k such that j + k < m. It follows 
that all distributions of order at most m and support x' x x" are linear combinations of the 
given tensor distributions. 

Our main tool is an adaptation of Theorem l6.2l 

Theorem 7.8. Given a test function <p and a cell K = K' x K" with K' = F' ,(X') and 
K" — F",,(X ,r ), there is a test function ip such that if/ = <p on K and 

m n 

||(A') J '(A")VL * C(m, n) £ £(^)^V^)"||(A0'(A'')V|| £ - w + e (7.7) 
for all < j < m and < k < n. 

Proof. The method for cutting-off a smooth function on a cell has already been described 
in the proof of Theorem l7.6l Since we cut off first in one variable and then in the other, the 
estimates from the proof of Theorem l6.2l mav be applied directly. Suppose that we cut off 
in the first variable and then in the second. Taking ( 16.41 ) for the Laplacian (A') k in the first 
variable on a fixed slice U' x \y"} and substituting from the second lines of both of ( 16.21 i 
and ( 16.31 ). gives 

m 

l|(A')VL (W}) =S CW^&tW^w* + e 

1=0 

provided j < m. In this calculation we used that the harmonic part of a function (which 
was denoted H l -j in the proof of Theorem 16.21 ) is bounded by the L°° norm of the func- 
tion because of the maximum principle, and we extracted the scaling factor r' ,fi' , of the 
Laplacian on K' = FL(X') using the same argument as in ( 16.51 ). 

The same estimate is true with the same proof when if/ is replaced by (A")V ar, d by 
(A") k <p. We use this fact when we repeat the estimate in the second variable, because in 
this case we are cutting off the function that was modified at the first step. A little algebra 
then produces the desired estimate. □ 

Theorem 7.9. Let T be a distribution supported at (x' , x") e X' x X". Suppose that x' is 
such that either Theorem \6.6\ or Theorem \6.13\ ma\ be used to identify the distributions with 
support at x' , and make the same assumption for x" . Then T is a finite linear combination 
of tensor products T' X T" where T' is supported at xf and T" is supported at x". 

Proof. In light of the preceding discussion and Theorem l7.8l it suffices to show that if the 
given tensor distributions vanish on a test function (f> then the right side of ( 17.71 ) may be 
made less than 2eby taking K sufficiently small. The proof of this estimate is elementary: 
we simply go from (x',x") to (y',y") by using two Taylor-like expansions, one in each 
variable. 

Since T has compact support it also has finite order m. It then seems reasonable that 
each of the terms T x T" should be made up of a T of order k < m and a T" of order at 
most m-k. Unfortunately we cannot prove this in general because our scaling estimates are 
insufficiently refined, as was explained in Remark lpTTl This result is true if the distributions 
at x' and x" are such that none have scaling exactly equal to that of the Laplacian (meaning 
that if they are as in Theorem l6 . 5 1 then there is no y s equal to a power of rp, and if they are 
as in Corollary 16. 121 then there is no y s equal to a power of fi lv ). Given the limitations of 
our estimates we must instead allow the possibility that T is order k + 1 and T" is order 
m — k + 1 . 

Suppose then that T x T"4> = for all T of order up to k + 1 and T" of order up 
to m - k + 1. It follows that the differential (D") m ~ k vanishes on the one-variable smooth 
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function T'0(x', •)• The same reasoning as was used at the beginning of the proofs of 
Theorem |62] and Theorem l6 .131 shows that then T'(f> = o(r'^„fi'^,,) m ~ k on the set fx'} x K", 
so in particular at (x',y"). 

We now wish to repeat the argument to go from (x',y") to (y',y"). Instead of hav- 
ing vanishing distributions in the first variable at (x',y") we have only estimates on then- 
size, which we use to estimate the size of (D')'"<p. Recall from Definitions 16 . 3 1 and 16 . 1 01 
that the differential (D')'"<p for the second variable on the cell K' consists of a harmonic 
function with coefficients obtained using distributions of order at most m + 1, as well as 
G'(D')"' ~'A'0. where G' is the Green's operator for the cell K' . The harmonic function is it- 
self made up of pieces (one for each k < m) with scaling bounded by (rL-fiL,) ( or an equiv- 
alent quantity involving f3' w , ) and coefficients obtained using distributions in the first vari- 
able with order at most k + 1 . The estimate of the previous paragraph says that these coeffi- 
cients are o(r'^,fx'^,,) m ~ k , so each term of the harmonic functions is o{(r' w ,^[ v ,) k (r",,,ii" v ,,) m ~ k ) 
on K' X {y"}. A similar argument applies to G'(D')" 1 ^ 1 A'<p, because the G' produces an 
extra factor of r' w ,fi' w ,, and the harmonic piece of (D') m _1 A'0 that has scaling (rL//,,) is 
obtained via distributions of order at most k applied to A'0, each of which is a distribution 
of order k + 1 applied to 0. We may repeat this reasoning inductively across the terms of 
{D')"'(f> to obtain a bound of the form 

|(D')">| = oEfatoH. (7.8) 

k=0 

on the set K' x {/'}. Since we also know (from ( 16.101 ) and ( 16.251 )) that 

|0-(D')'0|= O (r; v ^) m 

on K' x {/'} we conclude that the estimate ( 17.8b is also true for <p itself. The point y" € K" 
was arbitrary, so we have 

IML-c*) = u^KoVxri (7.9) 

*:=() 

Our working thus far has shown that if T x T"<p = for all T of order up to k + 1 and 
T" of order up to m - k + 1, then ( 17.91 ) holds. However, this assumption obviously implies 
that V x r"((A')'(A")'0) = for all V of order up to k + 1 - I and T" of order up to 
m — k+1- Of l + i< m and < k < (m — i — t). Thus ( 17.9b improves to 

m-i-l 

||(A')'(A")>L W = o[ J) (r>^)V>^) m - ! '-'-* l (7.10) 
Substituting into ( 17.7b for the cutoff of (f> yields 

m n m-i-l 

||(A'y(A")7L ^ e + °E Yy^^ } ^-^ k Z iC&tiC-iC-r***} 

1=0 i=0 j=0 

The simplest way to complete the argument is to choose K' and K" such that r' w ,/j,' w , and 
r w"Kv" are comparable, at which point all terms in the sum are bounded. It follows that the 
sum term is o(l) so can be made less than e by requiring that K' and K" are also sufficiently 
small. □ 
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8. Hypoellipticity 



An important question in the analysis of PDE is to identify conditions under which a 
distributional solution of a PDE is actually a smooth function. In Euclidean space, an 
archetypal example is Weyl's proof that a weak solution of the Laplace equation is actually 
C°°. In order to study these questions one uses the notion of hypoellipticity, which we may 
now define in the setting of fractafolds based on p.c.f. fractals and their products. We will 
not settle any of the questions about hypoellipticity here, but simply suggest some natural 
problems for which the distribution theory we have introduced is the correct setting. 

We first define the singular support of a distribution, which intuitively corresponds to 
those points where the distribution is not locally smooth. 

Defintion 8.1. A distribution T is smooth on the open set Qi c Q if there is u e £(£2i) 



Using Lemma [331 for the case of a single p.c.f. fractal, or the analogous result derived 
from Theorem 17. 7 1 in the product setting, we see that if T is smooth on Cli and on O2 then 
it is smooth on Qi U Q.2, thus there is a maximal open set on which T is smooth. 

Defintion 8.2. For a distribution T, Let Q7- be the maximal open set on which T is smooth. 
The singular support of T is the set 



Let P be a polynomial of order k on R m , so P{g) = 2w<* fl *£* where k = K\ . . . K m is a 
multi-index, \k\ - £ Kj is its length, and i- K — f] ■ Consider the linear differential oper- 
ator P(A) = P(A\ , . . . , A m ) on a product F] Xj of p.c.f. self-similar fractals Xj with Lapla- 
cians A;. It is clear that for any distribution T we have SingSppt(P(A)r) c SingSppt(T), 
because when T is represented by u e £(Qr) then P(A)T is represented by P(A)u. By anal- 
ogy with the Euclidean case, we define a class of constant coefficient hypoelliptic linear 
differential operators. 

Defintion 8.3. P(A) is called hypoelliptic if SingSppt(P(A)r) = SingSppt(r) for all T e 
£>'(fi). 

Given the importance of hypoelliptic operators in the analysis of PDE on Euclidean 
spaces, it is natural to seek conditions that imply hypoellipticity of an operator on a p.c.f. 
fractal or on products of p.c.f. fractals. We expect that if P(A) is elliptic then it should be 
hypoelliptic; it also seems possible that the celebrated hypoellipticity criterion of Hormander 
ifTTl Section 11.1] might imply hypoellipticity in the fractal case, though we do not expect 
conditions of this type to be necessary because of examples like that motivating Conjec- 
turelBIJI 

Defintion 8.4. For a polynomial P(^) = Yj\K\<k a K^ K , the principal part of P is Po = 
Yt\K\=k a K% K ■ P is called elliptic if Po{%) + for £ + 0; equivalently P is elliptic if there 
is c> so |P (£)| > c\$ k . We call P(A) elliptic if P{f v ... ,f m ) is elliptic. 

Remark 8.5. The above definition is consistent with the usual one in the case that X is a 
Euclidean interval rather than a fractal set, but they do not coincide because we are dealing 
with a smaller class of operators. Specifically, for such X the Laplacian is d 2 /dx 2 , so our 
class of operators {P(A)} is smaller than the usual collection of constant coefficient linear 
partial differential operators P(d/dxi, . . ., d/dx m ). Similarly our class of elliptic operators 
is a strict subset of the classical one. 



such that 




singSppt(r) = sppt(r) \ n T 
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Conjecture 8.6. If X is a p.c.f. fractal and P(A) is an elliptic operator on the product 
space X m , then P(A) is hypoelliptic. 

In the case that m = 1, all operators P(A) are elliptic, and they can all be shown to be 
hypoelliptic. Indeed, by factoring the polynomial we can reduce to the case of the linear 
polynomial A + c for some complex constant c. The hypoellipticity of A + c is readily 
obtained from the fact that on small cells (A + c) has a resolvent kernel that is smooth away 
from the diagonal, as may be seen by representing the resolvent as an integral with respect 
to the heat kernel or by applying results from lfT2l . 

Conjecture 8.7. A sufficient condition for the hypoellipticity of P(A) is that D" P(g) / P(^) — > 
as £, — > oo for any partial derivative D K with \k\ > (compare to Theorems 11.1.1 
and 11.1.3 of tTQj. 

In the Euclidean setting the above condition is necessary as well as sufficient, but we 
do not expect this to be the case on fractals. In essence, the idea is that hypoellipticity 
of P(A) should depend only on whether the principal part Po(A) is hypoelliptic, and that 
this is equivalent (on the Fourier transform side) to estimates when inverting the algebraic 
equation 

Po(A u A m )u(Ai,. . . , A m ) = f(A\, . . . , A m ) 
for any choice of (A\, . . ., A m ) with each Aj an eigenvalue of Ay. Since all of these Aj are 

negative, ellipticity of P(A) says that |Po(^i, ■ • ■ , A m )\ > c\A\-\ \-A m \ for such (Ai, . . . , A m ), 

and this is sufficient to show the Fourier transform u has faster decay than /, so u should be 
as smooth or smoother than /. However the ellipticity condition should only be necessary 
if the points (A\, . . . , A m ) are dense in the positive orthant : gj > 0). In (6] it is shown 
that this is not the case for the Sierpinski Gasket fractal; specifically it is shown that in the 
case m = 2, the points (Ai , A2) omit an open neighborhood of a ray in the positive orthant. 
It follows that there is a > and b < such that aA\ + bAi is not elliptic but -a/b lies in 
the omitted neighborhood, so \aAi + bAi\ > c\A\ + /I2I whenever Aj is an eigenvalue of Aj. 
Following |61 we call operators of this type quasielliptic. Given that quasielliptic operators 
satisfy elliptic-type estimates on the spectrum, it seems likely that they will have similar 
smoothness properties to elliptic operators; Lp estimates for these operators may be found 
in recent work of Sikora 1501 . 

Deflntion 8.8. The operator P(A) is quasielliptic if there is c > such that |Po(£)l ^ c|£| 
for all £ € {{Ai, . . ., A m ) : Aj is an eigenvalue of A ; }. 

Conjecture 8.9. The quasielliptic operators of [611 are hypoelliptic. 
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